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Abstract 


At the end of the inflationary stage of the early universe, profuse particle 
production leads to the reheating of the universe. Such explosive particle 
production is due to parametric amplification of quantum fluctuations for the 
unbroken symmetry case (appropriate for chaotic inflation), or spinodal in¬ 
stabilities in the broken symmetry phase (which is the case in new inflation). 
This mechanism is non-perturbative and depends on the details of the particle 
physics models involved. A consistent study of this mechanism requires a de¬ 
tailed analysis and numerical treatment with an approximation scheme that 
ensures energy (covariant) conservation and a consistent non-perturbative im¬ 
plementation. 

We study the 0{N) symmetric vector model with a quartic self-interaction 
in the large N limit, Hartree and resummed one-loop approximations (with 
= 1) to address the non-perturbative issues. The non-equilibrium equa¬ 
tions of motions, their renormalization and the implementation of the ap¬ 
proximations are studied in arbitrary spatially flat FRW cosmologies. A full 
description, analytically and numerically is provided in Minkowski space-time 
to illustrate the fundamental phenomena in a simpler setting. 

We give analytic results for weak couplings and times short compared to 
the time at which the fluctuations become of the same order as the tree level 
terms, as well as numerical results including the full backreaction. In the case 
where the symmetry is unbroken, the analytical results agree spectacularly 
well with the numerical ones in their common domain of validity. In the 
broken symmetry case, interesting situations, corresponding to slow roll initial 
conditions from the unstable minimum at the origin, give rise to a new and 
unexpected phenomenon: the dynamical relaxation of the vacuum energy. 


*Laboratoire Associe au CNRS UA280. 
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That is, particles are abundantly produced at the expense of the quantum 
vacuum energy while the zero mode comes back to almost its initial value. 
We obtain analytically and numerically the equation of state which in both 
cases can be written in terms of an effective polytropic index that interpolates 
between vacuum and radiation-like domination. 

The self-consistent methods presented in these lectures are the only ap¬ 
proaches, so far, that lead to reliable quantitative results on the reheating 
mechanism in the inflationary universe. These approaches take into account 
the non-linear interaction between the quantum modes and exactly conserve 
energy (covariantly). Simplified analysis that do not include the full backre- 
action and do not conserve energy, result in unbound particle production and 
lead to quantitatively erroneous results. 

For spontaneously broken theories the issue of whether the symmetry may 
be restored or not by the quantum fluctuations is analyzed. The precise crite¬ 
rion for symmetry restoration is presented. The field dynamics is symmetric 
when the energy density in the initial state is larger than the top of the tree 
level potential. When the initial energy density is below the top of the tree 
level potential, the symmetry is broken. 

Finally, we provide estimates of the reheating temperature as well as a 
discussion of the inconsistency of a kinetic approach to thermalization when 
a non-perturbatively large number of particles is created. 


I. INTRODUCTION 

Research activity on inflationary cosmologies has continued steadily since the concept of 
inflationary cosmology was first proposed in 1981 [|I|. 

It was recognized that in order to merge an inflationary scenario with standard Big 
Bang cosmology a mechanism to reheat the universe was needed. Such a mechanism must 
be present in any inflationary model to raise the temperature of the Universe at the end of 
inflation, thus the problem of reheating acquired further importance deserving more careful 
investigation. The original version of reheating envisaged that during the last stages of 
inflation when the universe expansion slows down, the energy stored in the oscillations of 
the infiaton zero mode transforms into particles via single particle decay. Such particle 
production reheats the universe whose temperature was redshifted to almost zero during the 
inflationary expansion 0]. 

It was realized recently that the elementary theory of reheating does not 

describe accurately the quantum dynamics of the fields when the oscillations of the infiaton 
field (zero mode) have large amplitude. 

Our programme on non-equilibrium dynamics in quantum field theory, started in 1992 
0, is naturally poised to provide a framework to study these problems. The larger goal of the 
program is to study the dynamics of non-equilibrium processes, such as phase transitions, 
from a fundamental field-theoretical description, by obtaining and solving the dynamical 
equations of motion for expectation values and correlation functions of the underlying four 


2 


dimensional quantum field theory for physically relevant problems: phase transitions and 
particle production out of equilibrium, symmetry breaking and dissipative processes. 

The focus of our work is to describe the quantum held dynamics when the energy den¬ 
sity is large. That is, a large number of particles per volume m~^, where m is the typical 
mass scale in the theory. Usual S-matrix calculations apply in the opposite limit of low 
energy density and since they only provide information on in —>■ out matrix elements, are 
unsuitable for calculations of time dependent expectation values. Our methods were natu¬ 
rally applied to different physical problems like pion condensates supercooled phase 

transitions PPP] , inhationary cosmology PJ5|J5|,P!3HT7| , the hadronization stage of the quark- 
gluon plasma |]^ as well as trying to understand out of equilibrium particle production in 
strong electromagnetic helds and in heavy ion collisions P,P,[r^ . 

When a large energy density is concentrated in one or few modes, for example the 
inhaton zero mode in inhationary cosmology, under time evolution this energy density will 
be transferred to other modes driving a large amplihcation of quantum huctuations. This, in 
turn, gives rise to profuse particle production for bosonic helds, creating quanta in a highly 
non-equilibrium distribution, radically changing the standard picture of reheating the post¬ 
inhationary universe PJT^. Fermionic helds are not very efficient for this mechanism of 


energy “cascading” because of Pauli blocking . 

The detail of the processes giving rise to preheating can be diherent depending on the 
potential for the scalar held and couplings to other helds involved, as well as the initial con¬ 
ditions. For example, in new inhationary scenarios, where the expectation value of the zero 
mode of the inhaton held evolves down the hat portion of a potential admitting spontaneous 
symmetry breaking, particle production occurs due to the existence of unstable held modes 
whose amplitude is amplihed until the zero mode leaves the instability region. These are 
the instabilities that give rise to spinodal decomposition and phase separation. In contrast, 
if we start with chaotic initial conditions, so that the held has large initial amplitude, par¬ 
ticles are created from the parametric amplihcation of the quantum huctuations due to the 
oscillations of the zero mode and the transfer of energy to higher modes. 

In these lectures we analyze the details of this so-called preheating process both an¬ 
alytically as well as numerically. Preheating is a non-perturbative process, with typically 
1/A particles being produced, where A is the self coupling of the held. Due to this fact, any 
attempts at analyzing the detailed dynamics of preheating must also be non-perturbative in 
nature. This leads us to consider the 0{N) vector model in the large N limit. This is a non- 
perturbative approximation that has many important features that justify its use: unlike 
the Hartree or mean-held approximation |P], it can be systematically improved in the 1/N 
expansion. It conserves energy, satishes the Ward identities of the underlying symmetry, 
and again unlike the Hartree approximation it predicts the correct order of the transition in 
equilibrium. 

This approximation has also been used in other non-equilibrium contexts . 

Our main results can summarized as follows . 


We provide consistent non-perturbative analytic estimates of the non-equilibrium pro¬ 
cesses occurring during the preheating stage taking into account the exact evolution of 
the inhaton zero mode for large amplitudes when the quantum back-reaction due to the 
produced particles is negligible i.e. at early and intermediate times. We also compute the 
momentum distribution of the number of particles created, as well as the ehective equa- 
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tion of state during this stage. Explicit expressions for the growth of quantum fluctuations, 
the preheating time scale, and the effective (time dependent) polytropic index dehning the 
equation of state are given in sec. IV and V. 

We go beyond the early/intermediate time regime and evolve the equations of motion 
numerically, taking into account back-reaction effects. (That is, the non-linear quantum held 
interaction). These results conhrm the analytic estimates in their domain of validity and 
show how, when back-reaction effects are large enough to compete with tree level effects, 
dissipational effects arise in the zero mode. Energy conservation is guaranteed in the full 
backreaction problem, leading to the eventual shut-off of particle production. This is an 
important ingredient in the dynamics that determines the relevant time scales. 

We also hnd a novel dynamical relaxation of the vacuum energy in this regime when 
the theory is in the broken phase. Namely, particles are produced at the expense of the 
quantum vacuum energy while the zero mode contributes very little. We hnd a radiation 
type equation of state for late times (p ~ | e) despite the lack of local thermodynamic 
equilibrium. 

Finally, we provide an estimate of the reheating temperature under clearly specihed (and 
physically reasonable assumptions) in a class of models. We comment on when the kinetic 
approach to thermalization and equilibration is applicable. 

iH 


There have been a number of papers (see refs. 


2^) dedicated to the analysis of 


the preheating process where particle production and back-reaction are estimated in different 
approximations |B3. Our analysis differs from other works in many important aspects. We 


emphasize the need of a non-perturbative, self-consistent treatment that includes backreac¬ 
tion and guarantees energy conservation (covariant conservation in the expanding universe) 
and the conservation of all of the important symmetries. Although analytic simplihed argu¬ 
ments may provide a qualitative picture of the phenomena involved, a quantitative statement 
requires a detailed numerical study in a consistent manner. Only a self-consistent, energy 
conserving scheme that includes backreaction effects can capture the corresponding time 
scales. Otherwise inhnite particle production may result from uncontrolled approximations. 

The layout of these lectures is as follows. Section II presents the model, the evolution 
equations, the renormalization of the equations of motion and introduces the relevant dehni- 
tions of particle number, energy and pressure and the details of their renormalization. The 
unbroken and broken symmetry cases are presented in detail and the differences in their 
treatment are clearly explained. 

In sections III through V we present a detailed analytic and numerical treatment of 
both the unbroken and broken symmetry phases emphasizing the description of particle 
production, energy, pressure and the equation of state. In the broken symmetry case, when 
the inflaton zero mode begins very close to the top of the potential, we hnd that there is a 
novel phenomenon of relaxation of the vacuum energy that explicitly accounts for profuse 
particle production through the spinodal instabilities and energy conservation. We discuss 
in section VI why the phenomenon of symmetry restoration at preheating, discussed by 
various authors |7| J2(]|J3^JS3[| is not seen to occur in the cases treated by us in ref. [&0] and 
relevant for new inhationary scenarios . 


A precise criterion for symmetry restoration is given. The symmetry is broken or un¬ 
broken depending on the value of the initial energy density of the state. When the energy 
density in the initial state is larger than the top of the tree level potential then the sym- 
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metry is restored When it is smaller than the top of the tree level potential, then it 
is broken and Goldstone bosons appear In the hrst case, the amplitude of the zero 


mode is such that V{rjQ) > ld(0) (all energy is initially on the zero mode). In this case the 
dynamics is very similar to the unbroken symmetry case, the amplitude of the zero mode 
will damp out, transferring energy to the quantum fluctuations via parametric amplihcation, 
but asymptotically oscillating around zero with a fairly large amplitude. 

In section VII we briefly discuss the amplitude expansion (linearizing in the held am¬ 
plitude) and compare with a full non-linear treatment in a model for reheating where the 
inhaton decays into a lighter scalar held TO . 


In section VIII we provide estimates, under suitably specihed assumptions, of the reheat¬ 
ing temperature in the 0{N) model as well as other models in which the inhaton couples to 
lighter scalars. In this section we argue that thermalization cannot be studied with a kinetic 
approach because of the non-perturbatively large occupation number of long-wavelength 
modes. 

Finally, we summarize our results and discuss future avenues of study in the conclusions. 


II. NON EQUILIBRIUM SCALAR FIELD DYNAMICS AT LARGE ENERGY 

DENSITIES 


Two essential parameters characterize the dynamics of quantum helds: the strength of 
the coupling A and the energy density in units of the typical mass m. If initially most of the 
energy is stored in one (or few) modes, the energy density is controlled by the amplitude 
of the expectation value of such mode(s) A = '/X <F/m. Usual held theory treatments 
consider the small amplitude limit A << 1 in which case the dynamics essentially reduces 
to the calculation of S'-matrix elements. The F-matrix describes the interaction of typically 
few particles in inhnite space-time. This is within the small amplitude limit even for high 
energies. 

We shall be concerned here with the non-perturbative regime in A = a/A $/m ~ 
(9(cxo). The crucial point is that non-linear ehects appear in such regime even for very small 

A. 

The small amplitude limit is also instructive to study as an initial condition prob¬ 

lem. In such regime, the held evolution equations linearize and can be solved explicitly by 
Laplace transform. Moreover, their solution can be interpreted using the S'-matrix language: 
contributions from particle poles, production thresholds for many-particle cuts, and so on. 

We consider the 0{N) vector model with quartic interaction in a cosmological spacetime 
with metric 

= dt^ — a^(t) 

Here, a{t) is the scale factor and t is the cosmic time coordinate. 

The action and Lagrangian density are given by. 


S = 


f 


d'^x £, 


C 


a^{t) 


1 -2 
-<F (x) 
2 ^ ^ 


1 ( V $( a ;))2 

2 a^it) 


U($(a;)) , 


5 











2 


( 2 . 1 ) 


y($) 


m 



2Nm?{t) \ 


N rn^ii) 

2A 


rn?‘{t) = rn? + ^ 7l{t) 


Here, TZit) stands for the scalar curvature. 


7^(^) = 6 


d{t) a?{t) 
a{t) a^{t) 


2 

where we have included the coupling ^ of <h(a:) to the scalar curvature since it will arise as 
a consequence of renormalization. The canonical momentum conjugate to <h(x) is, 


n(a;) = a?{t) <h(x). 


and the time dependent Hamiltonian is given by, 

H(t) = I A jlM + f + a^t) . 

In general, the system is in a mixed state described by a density matrix /5(<I)(.),$(.), t) in 
the Fock space. Here <I>(.) and <h(.) label the row and columns of the density matrix. Its 
time evolution is dehned by the quantum Liouville equation 

it^ = lH(t),p] ( 2 . 2 ) 

and we normalize it according to 

Tr p = 1 . 

The expectation value of any physical magnitude A is given as usual by 


< A >= Tr[p^] . 

The time evolution of all physical magnitudes is unitary as we see from eq. (|2.2|) . This 
implies that Von Neuman’s entropy 


S = Tr[p logp] . 


is conserved in time. 

In the present lectures we will restrict ourselves to translationally invariant situations. 
Namely, the order parameter 

< <F(x, t) > 

will be independent of the spatial coordinates x. 

There are two approximation schemes that have been used to study the non-equilibrium 
dynamics during phase transitions, each with its own advantages and disadvantages. The 

Igg,0,0,0 has the advantage that it can treat the dynamics of 
a scalar order parameter with discrete symmetry, while its disadvantage is that it is difficult 
to implement consistently beyond the lowest (mean held) level. The advantage of the large 


Hartree factorization 23 
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N approximation is that it allows a consistent expansion in a small 

parameter (1/A^) and correctly treats continuous symmetries in the sense that it implements 
Goldstone’s theorem.Moreover, the Hartree approximation becomes the resummed one-loop 
approximation for small values of A. Therefore, it may be a reliable approximation for the 
typical values of A in inflationary models. It should be noted that for spontaneous symmetry 
breaking, the large N limit always produces massless Goldstone bosons. 

Both methods implement a resummation of a select set of diagrams to all orders and 
lead to a system of equations that is energy conserving in Minkowski space time, and as will 
be shown below, satishes covariant conservation of the energy momentum tensor in FRW 
cosmologies. Furthermore, both methods are renormalizable and numerically implementable. 
Given that both methods have advantages and disadvantages and that choosing a particular 
scheme will undoubtedly lead to criticism and questions about their reliability, we use both, 
comparing the results to obtain universal features of the dynamics. 

In this section we introduce the 0{N) vector model, obtain the non-equilibrium evolution 
equations both in the large N and Hartree approximations, the energy momentum tensor 
and analyze the issue of renormalization. We will then be poised to present the analytical 
and numerical solutions as well as the analysis of the physics in the later sections. 

We choose the coupling A hxed in the large N limit. The held <I> is an 0{N) vector, 
<h = {a, tt) and if represents the iV — 1 “pions”. In what follows, we will consider two 
different cases of the potential (|2.1|) V{a,7r), with (m^ < 0) or without (m^ > 0) symmetry 
breaking. 

We can decompose the held a into its zero mode and huctuations x(T, t) about the zero 
mode: 

a{x,t) = ao{t) -f x(T,t) . 


The generating functional of real time non-equilibrium Green’s functions can be written 
in terms of a path integral along a complex contour in time, corresponding to forward and 
backward time evolution and at hnite temperature a branch down the imaginary time axis. 
This requires doubling the number of helds which now carry a label ± corresponding to 
forward (-I-), and backward (—) time evolution. The reader is referred to the literature 
for more details This generating functional along the complex contour requires the 

Lagrangian density along the contour, which is given by |p] 


>C[ao + x+,7f+]-£[cro + x , vf ] = |£[ao, 7f+] + 




2^^ ’ 2 a{ty 


, l(V7f )■ 





The tadpole condition {x'^{x,t)) 
P] and references therein. 


0 will lead to the equations of motion as discussed in 
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A. The Large N limit 


A consistent and elegant version of the large N limit for non-equilibrium problems can 
be obtained by introducing an auxiliary field (see for example This formulation has 

the advantage that it can incorporate the 0{1/N) corrections in a systematic fashion. Al¬ 
ternatively, the large N limit can be implemented via a Hartree-like factorization |]8|,[l5[| in 


which i) there are no cross correlations between the pions and sigma held and ii) the two 
point correlation functions of the pion held are diagonal in the 0(A —1) space of the remain¬ 
ing unbroken symmetry group. To leading order in large N both methods are completely 
equivalent and for simplicity of presentation we chose the factorization method. 

The factorization of the non-linear terms in the Lagrangian is (again for both ± compo¬ 
nents) : 

6 (x^) + constant 

^ 3 (x^) X 

(vf ■ tt)^ —^ 2 (tt^) — (vf^)^ -|- 0{1/N) 

vfV-(^^)x^ + vf^(x^) 

^^x ^ (^^)x 


To obtain a large N limit, we dehne 


N-l 


(2.3) 


i{x,t) = ■ ■ ,1) ; a„{t) = 4>{t)VN 

where the large N limit is implemented by the requirement that 

^ Oil) , (x^) ^ 0(1) , 0 ^ 0(1). 

The leading contribution is obtained by neglecting the 0(1/A) terms in the formal large N 
limit. 


8IL 

>C[cro x"^, 7f+] - £[ao + X“, 7f“] = <( £[ao, 7f+] -h ^X"^ + a^(^) 


1._ 1 (v^ 

2^^ ^ 2 a(f)2 


-(# 

2^ 


~ \ ^I(V^ ~ (^^"(^o,7f+)x+^ + ^V"[3](cro,ir+)(x+)^ + 7f+)(x+)^) 


The resulting Lagrangian density is quadratic, with a linear term in x : 

1 (Vx^)" , 


>C[cro x"^, - >C[cro X , ] = <^ a^(t) 


l(y+)2- 

2^^ ’ 2 a(f)2 


.vr — 


2 a(t)2 
1 


- x^v'it) 


^Mlit)ix^)^ - ^-M|(t)(vf+)2} - {(x+ ^ X“) , ^ vf-)} 


(2.4) 



















where, 


V'm),t) = s/Nm 




Mlit) = m{tf + 

where m{tY is dehned in eq.( p.l|) . Note that we have used spatial translational invariance 
as behts a spatially flat FRW cosmology, to write 

When the initial state is in local thermodynamic equilibrium at temperature Ti, the flnite 
temperature non-equilibrium Green’s functions are obtained from the following ingredients 

Glit.t') = l^{fk{t)fl(t')[l + nk]+nk fk{G)fk{t)} 

^ {fk{t')fk{t)[l + nk]+nk fk{t)fk{t')} 


Efc 


where Uk = {e'^i — 1) 

The Heisenberg held operator 'ip{x,t) can be written as 




Sk 1 

(27r)3 ^2 L 


H fki't) + 4 fkit) e 


t f*i 


-ik-x 


( 2 . 6 ) 


where ak , al are the canonical destruction and annihilation operators. 

The evolution equations for the expectation value 0(f) and the mode functions fk{t) can 
be obtained by using the tadpole method |p[ and are given by: 


0(f) + 3if0(f) + m(f)^0(f) + ^0^(^) + ^ 0(^) (0^(^))i3 = 0 , 


with the mode functions. 


( 2 . 6 ) 


^ + STT—a;fc(f) /^(f) — 0, 


(2.7) 


and the effective frequencies. 


“ ~TF7\ + ) 


a2(f) 


where the effective mass takes the form. 


M2(f) = m{tf + ^0^(f) + ^(0^(f)) 


B ■ 


( 2 , 8 ) 


Here, the bare quantum fluctuations are given in terms of the mode functions by , 
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(2.9) 




d^k |/fc(t)p 

(27r)3 2 


coth 


.2Ti. ■ 


At this stage we must provide the initial conditions on the mode functions fk{t)- As 
mentioned above our choice of initial conditions on the density matrix is that of local ther¬ 
modynamic equilibrium for the instantaneous modes of the time dependent Hamiltonian at 
the initial time. Therefore we choose the initial conditions on the mode functions to repre¬ 
sent positive energy particle states of the instantaneous Hamiltonian at f = 0, which is the 
initial time. Therefore our choice of boundary conditions at f = 0,is 

fk{0) = ; MO) = -^]/m ; Wk = 

where the mass Mq determines the frequencies a;fc( 0 ) and will be obtained explicitly later. 
With these boundary conditions, the mode functions /^(O) correspond to positive frequency 
modes (particles) of the instantaneous quadratic Hamiltonian for oscillators of mass Mq. 
The initial density matrix, at time f = 0 is thus chosen to be that of local thermodynamic 
equilibrium at the temperature Tj for these harmonic modes. The fluctuations x(a;, t) obey an 
independent equation, that does not enter in the dynamics of the evolution of the expectation 
value or the vf fields to this order and decouples in the leading order in the large N limit [§ . 

It is clear from the above equations that the Ward identities of Goldstone’s theorem are 
fulhlled. Because V'{(p{t),t) = y/N(p{t)MM)^ whenever V'{(p{t),t) vanishes for 0 7 ^ 0 then 
= 0 and the “pions” are the Goldstone bosons. This observation will be important in 
the discussions of symmetry breaking in a later section. 

Since in this approximation, the dynamics for the vf and x fields decouple, and the 
dynamics of y does not influence that of 0 , the mode functions or (' 0 ^), we will only con¬ 
centrate on the solution for the vf helds. We note however, that if the dynamics is such that 
the asymptotic value of 0 7 ^ 0 the masses for y and the “pion” multiplet vf are different, and 
the original 0{N) symmetry is broken down to the 0{N — 1) subgroup. 


B. The Hartree and the One-loop Approximations 

To implement the Hartree approximation, we set iV = 1 and write, 

with, 

0(f) = ($(T,f)) ; (0(f,f)) = 0, 

where the expectation value is dehned by the non-equilibrium density matrix specihed below, 
and we have assumed spatial translational invariance, compatible with a spatially flat metric. 
The Hartree approximation is obtained after the factorization, 

3 0(T,f), 


M{x,t) 6 {M{x,t)) 0^(f,f) -3 {M{x,t))^, 
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where by translational invariance, the expectation values only depend on time. In this 
approximation, the Hamiltonian becomes quadratic at the expense of a self-consistent con¬ 
dition. 

At this stage we must specify the non-equilibrium state in which we compute the expec¬ 
tation values above. In non-equilibrium held theory, the important ingredient is the time 
evolution of the density matrix p{t) (see and references therein). This density matrix 
obeys the quantum Liouville equation (|2.2|) whose solution only requires an initial condition 
p{U) [|2],|li|,|,|n|,|T5|,i. The choice of initial conditions for this density matrix is an issue 
that pervades any calculation in cosmology. Since we want to study the dynamics of the 
phase transition, it is natural to consider initial conditions that describe the instantaneous 
modes of the time dependent Hamiltonian as being initially in local thermodynamic equi¬ 
librium at some temperature Tj > Tc. Given this initial density matrix, we then evolve it 
in time using the time dependent Hamiltonian as in ^ or alternatively using the complex 
time path integral method as described in P^P^ , P!^ J8|J9| JTI|JT3| J^ . 

Following the steps of references P,| 


we hnd the equation of motion for the 


expectation value of the inflaton held to be, 

\ Q\ 

0(t) + 3H^{t) + M^(j){t) + + Y (V’^(^))b = 0 


( 2 . 10 ) 


The bare quantum huctuations {'ijj‘^{t))B are obtained from the coincidence limit of the non¬ 
equilibrium Green’s functions, which are obtained from the mode functions obeying. 


d 2 , . 


Mt) = 0 , 


( 2 . 11 ) 


with the effective frequencies. 


where 


^lit) — + M^{t) , 


= mity + + y (V’“(l)> 


( 2 . 12 ) 


Notice the only difference between Hartree and large N limits: a factor 3 in front of 
4>‘^{t) -f {'ijj‘^{t)) in the effective mass squared for the mode functions as compared to the 
equation for the zero mode. In particular, when 0(t) = 0 corresponding to a phase transition 
in absence of biased initial conditions, both descriptions yield the same results (up to a trivial 
rescaling of the coupling constant by a factor 3). 

The equal time correlation function is given in terms of the mode functions as 0J3,n, 


mt)) = 


|/fc(t)p 

(27r)3 2 


coth 


2T,; 


(2.13) 


The initial conditions are chosen to rehect the same physical situation as in the large 
N case, that is, the instantaneous particle states of the Hamiltonian at t = 0 are in local 
thermodynamic equilibrium at some initial temperature higher than the critical value. Thus, 
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as in the large N case but with modified frequencies, the initial conditions at t = 0 are chosen 
to describe the instantaneous positive energy states, 

/,(0) = ^ ; /fc(0) = -l^/Wk ; Wk = . (2.14) 

We have maintained the same names for the mode functions and Mq to avoid cluttering 
of notation; their meaning for each case should be clear from the context. Notice that 
the difference between the Hartree and large N case is rather minor. The most significant 
difference is that, in the equations for the zero modes, the Hartree case displays a factor 3 
difference between the tree level non-linear term and the contribution from the fluctuation 
as compared to the corresponding terms in the large N case. The equations for the mode 
functions are the same upon a rescaling of the coupling constant by a factor 3. 


A re-summed one-loop approximation is obtained by keeping only the leading quantum 
corrections. That is, the first non-trivial contribution in A. Such approximation can be 
worked out by taking the expectation value of the evolution equations of the field operator 
<h to first order in A. At this stage, we can straightforwardly obtain the resummed one- 
loop evolution equations from the Hartree equations for small A. Just notice that 
is multiplied by A in the zero mode equation (|2.10|) . Therefore, to leading order in A, we 


can neglect the term | A (V’^(t)) in the mode equations (p.ll|) . In summary, the resummed 
one-loop evolution equations take the form 

A 3A 

m + 3iJ0(f) -f + -03(f) + - 0(f) (02(f))s = 0 , 


— 3H— 

df^ df a^(f) 


^ +m(f)2 + y02(f) 


fk{t) = 0 . 


The bare one-loop quantum fluctuations are obtained by inserting the one-loop 


modes fk{t) into eq. (|2.13|) . 

It must be stressed, however, that a numerical implementation of the set of equations 
above, represents a non-perturbative treatment, in the sense that the (numerical) solution will 
incorporate arbitrary powers of A. A naive perturbative expansion in A is bound to break 
down due to secular terms whenever resonances are present as is the case in parametric 
amplification. A resummation of these secular terms as obtained via a numerical integration 
for example corresponds to a non-trivial resummation of the perturbative series. 

These resummed one-loop equations are slightly simpler than the large N or Hartree 
equations. For small values of A as in inflationary models (where A ~ 10“^^) the resummed 
one-loop approximation provides reliable results TO . 


C. Renormalization in Cosmological Spacetimes 


We briefly review the most relevant features of the renormalization program in the large 
N limit that will be used frequently in our analysis. The Hartree case follows upon trivial 
changes. For more details the reader is referred to [|l4| , ^p! ^ 3- 
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In this approximation, the Lagrangian is quadratic, and there are no counterterms. This 
implies that the equations for the mode functions must be hnite. This requires that 

where the subscripts B, R refer to bare and renormalized quantities, respectively. Dehning 

ipk{t) = fk{t) , ipk{0) = -^= , (y9fe(0) = -i\[Wk 

V trfc 


(with a(0) = 1). 

The functions (pk{t) satisfy the Schrodinger-like differential equation 


(R 3 / a 1 \ P , .r2 / \ \ 

— - x(- + x-l+^ + M\t) Mt) = 0 


dR 2 \ a 2 a? 


a\t) 


In order to derive the large k behaviour, it is convenient to write the ^k{t) as linear combi¬ 
nations of WKB solutions of the form 

(Pfc(^) = R-k exp [ Rk{t')dt' + Bk exp [ Rl.{t')dt' 

Jo Jo 

with Rk(t) obeying a Riccati equation P] and the coefficients Ak , Bk are fixed by the initial 
conditions. After some algebra we find , 


m)\^ = 


1 1 
+ 


ka^{t) 2k‘^a?{t) - 
1 


H‘^R))-B{t) 


d r 


8a{tR k^ 


B{t)[3B{t)-2HR0)]+a{t)— a{t)B{t) + Do \ + 0{l/kR 


\m\^ = 


k 


+ 


a‘^{t) ka^{t) 
1 


H\t) + 


hf 2 (o) 1 

2a?(t) ^ 2 


dt L 


M^{t) 


m' 


— ait^Bit) + 3a{t)d(t)B{t) — 4a^(t)i?(t) 


8a{tR k? 

+ 2H\Q)[2d^{t) + Bit)] + Do} + 0(1/P 

where we defined BR) as 

R(t) = a^(t) [ M^(t) 


(2.15) 


m ' 

6 


in terms of the effective mass term for the large N limit given by (|2.8|) and the Hartree case, 
eq. (|2.12|) . The constant Dq depends on the initial conditions and is unimportant for our 
analysis. 

Using this asymptotic forms, we obtain [R|,O,&0,IT^ the following renormalized quantities 


m|(t) + 


Xb XR Xb /A\ a? {to) 


+ 


IQir'^ a‘^{t) Ibvr^ \kJ a?it) 


In — 


= m 


R 


1 + 


A 


B 


IGtt^ 


Inf- 

K 
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(2.16) 


Xb — 


= ^R + 


X 


R 


Af 


167r^ 

IGtt^ 




In 


A 


R 


^R 


m?)R = 
e{k 


d?k 



rhFM 

1 2 

mil 


R] 


AkX a^{t) 


(27r)3 

-H‘^{Q) + a^{t) ( M2(t) 


2k a^{t) 


m ' 

6 


We have introduced the (arbitrary) renormalization scale r. The conformal coupling .^ = 1/6 
is a fixed point under renormalization In dimensional regularization the terms involving 
are absent and In A is replaced by a simple pole at the physical dimension. Even in such 
a regularization scheme, however, a time dependent bare mass is needed. The presence of 
this new renormalization allows us to introduce a new renormalized mass term of the form 

Q 

a?{t) 

This counterterm may be interpreted as a squared mass red-shifted by the expansion of the 
universe. However, we shall set p = 0 for simplicity. 

At this point it is convenient to absorb a further finite renormalization in the definition 
of the mass and introduce the following quantities: 


A 


u 

T = \MR\t , q = 


\Mr\ 


n =-^ 

\Mr\ 


T = 


\Mr\ 


7]\t) = 

g^ir) = 


X 


R 


2\Mr\^ 

Xr 

2\Mr[^ 




mt))R-mo))R 


S(0) = 0 


9 = 


X 


R 


(2.17) 

(2.18) 


87r2 ’ 

^qi^) = ^/\MR\ fk{t) ■ 

For simplicity in our numerical calculations later, we will chose the renormalization scale 
R = \Mr\. The evolution equations are now written in terms of these dimensionless variables, 
in which dots now stand for derivatives with respect to r. 


III. THE RENORMALIZED EVOLUTION EQUATIONS AND THE 
ENERGY-MOMENTUM TENSOR 

From now on we focus our analysis on the case of Minkowski space-time with the aim of 
understanding the fundamental phenomena in a simpler setting. The case of cosmological 
spacetimes is presented in refs. B0. 

Let us summarize here the renormalized field equations in the Hartree and large N 
approximations that will be solved numerically and analytically in Minkowski spacetime. 
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A. Unbroken Symmetry 


In this case M|. = and in terms of the dimensionless variables introduced above 

we hnd the following equations of motion: 


g + g + g^ + g g^r) S(r) = 0 


(3.1) 

^ + 1 + girf + g S(r) 

^q{r) = 0 , 

(3.2) 

9^g(0) — r— , V^g(O) — i \l~^q 

W i ig 

(3.3) 

^(0) = go , g{0) = 0 


(3.4) 


Hence, Al^(r) = 1 + + g S(r) plays the role of a (time dependent) renormalized 

effective mass squared. 

As mentioned above, the choice of hlg determines the initial state. We will choose these 
such that at f = 0 the quantum fluctuations are in the ground state of the oscillators at the 
initial time. Recalling that by definition (y'S(O) = 0, we choose the dimensionless frequencies 
to be 


(f + ^ + 'll- 

The Wronskian of two solutions of (|3.2|) is given by 

[^q, ^q] = 2 * , 


while 5 'S(r) is given by 


POO ( 


^q 

2T 


2g3 


-ril + 7f{T)+gT.{T) 



(3.5) 


(3.6) 


B. Broken Symmetry 

In the case of broken symmetry = — |M^| and the field equations in the N = oo 
limit become: 


g -g + T]^ + g g{T) E(r) = 0 
^ + q‘ -1+ t)(r)^ + g E(r) 


= 0 


(3.7) 

(3.8) 


where S(r) is given in terms of the mode functions (pqir) by the same expression of the 
previous case, ( |3.6| ). Here, = —1+ g{T)‘^+g S(r) plays the role of a (time dependent) 

renormalized effective mass squared. 
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The choice of boundary conditions is more subtle for broken symmetry. The situation of 
interest is when 0 < < 1, corresponding to the situation where the expectation value rolls 

down the potential hill from the origin. The modes with <1 — rj^ are unstable and thus 
do not represent simple harmonic oscillator quantum states. Therefore one must chose a 
different set of boundary conditions for these modes. Our choice will be that corresponding 
to the ground state of an upright harmonic oscillator. This particular initial condition 
corresponds to a quench type of situation in which the initial state is evolved in time in an 
inverted parabolic potential (for early times f > 0). Thus we shall use the following initial 
conditions for the mode functions: 


V?g( 0 ) = 


0 g(O) = -i Jvtq 


+ ^ +Vo for <qI = ^-Vo 

-^ +Vo ioT > ql ; 0 < < 1 . 

along with the initial conditions for the zero mode given by eq. (|3.4|) . 


(3.9) 


(3.10) 


C. Particle Number 


Although the notion of particle number is ambiguous in a time dependent non-equilibrium 
situation, a suitable dehnition can be given with respect to some particular pointer state. 
We consider two particular dehnitions that are physically motivated and relevant as we will 
see later. The hrst is when we dehne particles with respect to the initial Fock vacuum state, 
while the second corresponds to dehning particles with respect to the adiabatic vacuum 
state. 

In the former case we write the spatial Fourier transform of the fluctuating held 'ijj{x,t) 
in (|2.3|) and (p.5|) and its canonical momentum n(T, t) as 


'ipkit) 

nfc(t) 


1 

V2 

1 

V2 


akfk{t) + aljkit) 

ttkfkit) + alJkit) 


with the time independent creation and annihilation operators, such that annihilates the 
initial Fock vacuum state. Using the initial conditions on the mode functions, the Heisenberg 
held operators are written as 


V’fc(i) = Id \t) -0^(0) U{t) 

Uk{t)=U-\t) Uk{0)U{t) 

dk{t) = U~\t) ttk U{t) 


1 

7K 


ak{t) + alfe(t) 




with U{t) the time evolution operator with the boundary condition W(0) = 1. The Heisen¬ 
berg operators dk{t) , d\{t) are related to ak, a], by a Bogoliubov (canonical) transformation 
(see reference [H for details). 
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The particle number with respect to the initial Fock vacuum state is dehned in term of 
the dimensionless variables introduced above as 


= (il(t)at(t)} 


= 2 ^.( 0 ) 




1 

+ 4 








1 

2 


(3.11) 


The initial occupation number iVg(O) exhibits a thermal distribution 

NM = ^ 


flq/r _ I ’ 


according to the initial temperature T. The particle number is expressed in eq. (|3.11|) as the 
sum of two contributions: the hrst term is the spontaneous production (proportional to the 
initial thermal occupation) and the second is the induced production (independent of it). 

It is the dehnition ( |3.11|) of particle number that will be used for the numerical study. 
In order to dehne the particle number with respect to the adiabatic vacuum state we 
note that the mode equations (|3.2| , |3.8|) are those of harmonic oscillators with time dependent 
squared frequencies 

± 1 + h^(T) + gS(r) 

with + for the unbroken symmetry case and — for the broken symmetry case, respectively. 
When the frequencies are real, the adiabatic modes can be introduced in the following 
manner: 


V’fc(t) 

nfc(t) 


2uk{t) 




where now akit) is a canonical operator that destroys the adiabatic vacuum state, and is 
related to ak, al by a Bogoliubov transformation. This expansion diagonalizes the instan¬ 
taneous Hamiltonian in terms of the canonical operators a{t) , a^{t). The adiabatic particle 
number is 


= (al(t)at(t)} 


= 2 


a;q(r) 


1 

+ 4 






(3.12) 

1 

' 2 ■ 


As mentioned above, the adiabatic particle number can only be dehned when the fre¬ 
quencies LOgir) are real. Thus, in the broken symmetry state they can only be dehned 

for wave-vectors larger than the maximum unstable wave-vector, k > ku = 

These adiabatic modes and the corresponding adiabatic particle number have been used 
previously within the non-equilibrium context [0 and will be very useful in the analysis 

Notice that 


of the energy below. Both dehnitions coincide at r = 0 because a;g(0) = klg. 

= A(j(0) = 0 if we choose zero initial temperature. (We considered a non-zero initial 


temperature in refs. [§,|1). 
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D. Energy and Pressure 


The energy-momentum tensor for this theory in Minkowski spacetime is given by 






(3.13) 


Since we consider translationally as well as rotationally invariant states, the expectation 
value of takes the fluid form 

Since we consider translationally and rotationally invariant states, the expectation value 
of the energy-momentum tensor takes the fluid form p =< P > /NV 

NV p{t) = < T^\x) >=< T^\x) >=< T^^{x) >=< ^ - T°°(a:) > , 

o 

with all off-diagonal components vanishing. 

Hence, 

p{t) + 0 > 

takes a particularly simple form. 

Using the large N factorization ( |2.3| - |2.3|) we find the energy density operator for zero 
initial temperature (T = 0) to be. 


E 1)0. X 1 0,0/x ^ id, ^ 1 f 

— = -0 (t) +-m 0 (f) +-0 W + 2 /^ 

— + linear terms in 0 -|- C>(l/iV) 

8 

ujlit) = e + m^ + ^0^(0 + ^(0^(t)). 




Analogous expressions can be derived for the energy in the Hartree approximation. 

The generalization to non-zero initial temperature is straightforward. 

Taking the expectation value in the initial state and the infinite volume limit (V —> cxd) 
and recalling that the tadpole condition requires that the expectation value of 0 vanishes, 
we hnd the expectation value of the bare energy to be 


Eharp. 


2\Mr\' 




0 ^ + girY + -pirY 


+ ^ dq I 00r) 


+ (g^ + 1 + virY) I (p0r) P]} + ^ {Y^Y)) 


2 

B ) 


(3.14) 


where is given by eq. (|2.9|) . It is easy to see that E^are = 0 using the bare equations 

of motion (|2.6| - |2.7|) . It is important to account for the last term when taking the time 
derivative because this term cancels a similar term in the time derivative of 

We want to emphasize that the full evolution of the zero mode plus the back-reaction 
with quantum fluctuations conserves energy (covariantly in expanding cosmologies). Such 
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is obviously not the case in treatments of reheating in the literature in which back-reaction 
effects on the zero mode are not taken into account in a self-consistent way. Without energy 
conservation, the quantum fluctuations grow without bound. In cosmological scenarios 
energy is not conserved but its time dependence is not arbitrary; in a hxed space-time 
background metric it is determined by the covariant conservation of the energy momentum 
tensor. There again only a full account of the quantum back-reaction will maintain covariant 
conservation of the energy momentum tensor. 

We hnd for the sum of bare energy plus pressure, 




bare 


Eharp — 


2\Mr[- 


A 


R 


^fi^ + gjq^dq I 0g(r) 1^I (Pg(r) p I . (3.15) 


It is clear that the integrals in eq. (p.l4|) and ( |3.15|) are divergent. 

In the previous section we have learned know how to renormalize ('0^(t))s, the renor¬ 
malized quantum fluctuations are denoted by S(r) [see eqs. ( |2.17|) and (|3.6[)]. 

In order to renormalize E^are and p{T)bare we need to use the large q behaviour of the 
mode functions Pqir) (|2.15|). In Minkowski spacetime this large q behaviour reduces to 




12 


mHt) 1 

2q^ 


12 9^°° 


= q + 


8 q^ 
1 


M^(t) _ 

2 q 8 q^ 


3 M‘(t) + 


+ o(«-') 

+ Oig-O) . 


(3.16) 


where M^(t) = ±1 -I- -|- g S(r) is the renormalized effective mass squared. 

We then subtract these asymptotic behaviours inside the integrand of eqs. (|3.14|) and 
( 3.15| ) in order to make the integral hnite. We hnd the following expression for the renor¬ 
malized energy setting A = cx) : 


E = 

■^rpT). 


Ai? 12 

2 I |2 


+ rjirf +-Viry 


+ U^q^dq [\pq{r) ^ 
2 JO L 


+ q I Pqi^) 

. 9 


-2q 


diq-y ,.4 


4g3 


M\t) 


E(r) 


1 + jEIt) 


It is easy to see that Eren is finite. Moreover, it is conserved. That is, we hnd that 
Eren = 0 using eqs. (|2.6|) and (p.7|) . 

We hnd that aside from the time independent divergence that is present in the energy 
the pressure needs an extra subtraction 


1 V 2 / 

6g3 dr^ 

compared with the energy. Such a term corresponds to an additive renormalization of the 
energy-momentum tensor of the form 

= A - d^d‘^)Syx) 
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with A a (divergent) constant |^. Performing the integrals with a spatial ultraviolet cutoff, 
and in terms of the renormalization scale k introduced before, we hnd 


A = -^lnl-] 

12 

In terms of dimensionless quantities and after subtracting a time independent quartic 
divergence, we finally find setting A = oo, for the renormalized energy plus pressure 


p(r) 


A'T’prj 


2\Mr\^ 


X 


R 


+ g 


dq 




- 3 ^- 


M^t) 9{q-K) 


3q 


12 q^ dr"^ f 


Ad^(r) 


In order to obtain a better insight on this quantum conserved energy it is convenient to 
write eq. (|3.14|) as 


Ebare 


2\MrY- 




?7^ + r7(r)2 + -p(r)' 

A 


+ ^ dq I (p,(r) 


+ I ‘dgM P]} - o > 


Then, we get using eq. (|3.12|) , 

[|v>,(t)|^+h)J(t)|v>,(t)|^] =irc; + 2^ 9^ w,(t) (y“'(r) + 1) , 


(3.17) 


1 rQu r 

"^= 2/0 <l^d<l[\^q{r)\^+u^^,{r)Mr)(^ 


where A is a spatial upper momentum cutoff, taken to inhnity after renormalization. In 
the broken symmetry case, eu is the contribution to the energy-momentum tensor from the 
unstable modes with negative squared frequencies, — pg] and is the 

adiabatic particle number given by eq. (|3.12|) . For the unbroken symmetry case eu = 0 and 

Qu = 0 . 

This representation is particularly useful in dealing with renormalization of the energy. 
Since the energy is conserved, a subtraction at r = 0 suffices to render it hnite in terms 
of the renormalized coupling and mass. Using energy conservation and the renormalization 
conditions in the large N limit, we hnd that the integral 


POO 

/ rfdqu,,(T)Nf(T) 


is hnite. This can also be seen from the asymptotic behaviors (B.lbl). We get from eqs. 
( gr^ ) and (|3T6D , 

NfU) ’3” 0(1) . 

All ultraviolet divergences are contained in the last term of eq.( p.l7|) . That is. 
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fA 

/ dq Ug{T) = 


1 


^ log(2A) + ^ M* 


(3.18) 


Qu\/ql + M^ (M^ + 2q^) - log (q^ + ^Jql + M^ 


In terms of dimensionless quantities, the renormalized energy density is, after taking 
A —> cxo: 


^ 2\Mn\Uv^ 1.,, 23 A.2/ 3 M\t) + 1 

<1 ^ + _(±i + 7j^)M\t) - Y -+ 9 


X 


R 


^F{r) + ^J^{vo)M^{t) 


qn 


{ql + M\t)) + ^ql + M'^{T)^ + y M‘^{T)^ql + 


M\t) M\t) 
32 ^ 8 


In 


qu + ^ql + M'^ir) 


+ C^M 


(3.19) 


where. 


Sf^t) 


1 

2 




q^dq |0/ + o^J(r)|(^/ +2 / q^dq ujq^r) Nf^r) 


' Qu 


±1 + t]‘^{t) + gJ:{T) , o^2(r) = + Af^(r) , 


(3.20) 

(3.21) 


Here the lower sign and = \Jl — apply to the broken symmetry case while the upper 
sign and g^ = 0 correspond to the unbroken symmetry case. The constant J^{'qo) is dehned 
as. 


POO 

J^iVo)= q^dq 
Jo 



2g3 


0 (g- 1) 


(3.22) 


The constant C^{riQ) is chosen such that Eren coincides with the classical energy for the zero 
mode at r = 0. The quantity M.{t) is identihed as the effective (dimensionless) mass for 
the “pions”. 

In the unbroken symmetry case (uper sign) we hnd 


J^{m) 


1 + ^g 

4 


1 + log 



and 


C 


\vo) = —^{i + vl) JAvo) 


We hnd using the renormalized eqs. ( ^rl) , (|3.2|) , (|3.6|) , (1^) and (|3.8|) , that the renor¬ 
malized energy Eren is indeed conserved both for unbroken and for broken symmetry. 


Let us now make contact with the effective potential which is a quantity dehned for time 
independent expectation value of the held. That is, for constant q. 

We recognize in eq. ( 3.19 ) that the sum of terms without Sp for g^ = 0 coincide with 
the ehective potential in this approximation. These arise from the ‘zero point’ energy of the 
oscillators in (|3.17|) . That is, for //(r) = r/oi 
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Veffim) 


2\Mr\" 


, (3.23) 


' M ^-1 


+ 9 






32 ■ 8 


+ 


Notice that A4^(t) = A4^ = ±1 + t^q for a time independent order parameter 7](t) = rjo as 
it follows from eq. (|3.6| ). 

In the broken symmetry case the term ejr describes the dynamics of the spinodal in¬ 
stabilities [Q since the mode functions will grow in time. Ignoring these instabilities and 
setting = 0 as is done in a calculation of the effective potential results in an imaginary 
part. In the unbroken symmetry (g„ = 0) case the sum of terms without Sjp give the effective 
potential in the large N limit. The term cannot be obtained in a purely static calcu¬ 
lation. Such term describes the profuse particle production via parametric amplihcation, 
the mode functions in the unstable bands give a contribution to this term that eventually 
becomes non-perturbatively large and comparable to the tree level terms as will be described 
in detail below. Clearly both in the broken and unbroken symmetry cases the effective po¬ 
tential misses all of the interesting dynamics, that is the exponential growth of quantum 
fluctuations and the ensuing particle production, either associated with unstable bands in 
the unbroken symmetry case or spinodal instabilities in the broken symmetry phase. 

The expression for the renormalized energy density given by ( |3.19[ - PT^TD differs from the 
effective potential in several fundamental aspects: i) it is always real as opposed to the 
effective potential that becomes complex in the spinodal region, ii) it accounts for particle 
production and time dependent phenomena. 

At this stage we can recognize why the effective potential is an irrelevant quantity to 
study the dynamics. 

The effective potential is a useless tool to study the dynamics precisely because it misses 
the profuse particle production associated with these dynamical, non-equilibrium and non- 
perturbative processes. 


IV. THE UNBROKEN SYMMETRY CASE 


The full resolution of the large N or Hartree equations needs a numerical treatment 
T^ . However an analytic treatment can be performed for early times while the non¬ 


linear effects are still small. 


A. Analytic Results for large N 


In this section we turn to the analytic treatment of equations ( |3.1|) , ( |3.2| ) and p.fj) in 
the unbroken symmetry case. Our approximations will only be valid in the weak coupling 
regime and for times small enough so that the quantum fluctuations, i.e. ^^(r) are not 
large compared to the “tree level” quantities. We will see that this encompasses the times 
in which most of the interesting physics occurs. 

Since H(0) = 0, the back-reaction term gH(r) is expected to be small for small g during 
an interval say 0 < r < ri. This time ri, to be determined below, determines the relevant 
time scale for preheating and will be called the preheating time. 


22 













During the interval of time in which the back-reaction term 5 'S(r) can be neglected 

fj + Tj + rf = Q . 


eg . (13.11) reduces to 


The solution of this equation with the initial conditions (p.4D can be written in terms of 
elliptic functions with the result: 


Vi.r) 


k 



Vo 

Y ^ 2(1 t]q) 


(4.1) 


where cn stands for the Jacobi cosine. Notice that ^^(r) has period 4lu = AKik)/ 
where K{k) is the complete elliptic integral of hrst kind. In addition we note that since 


7]{t + 2uj) = -7 ]{t) , 

if we neglect the back-reaction in the mode equations, the ‘potential’ (—1 —? 7 ^(r)) is periodic 
with period 2tn. Inserting this form for //(r) in eq. (|3.2|) and neglecting 5 'D(r) yields 

^ + 1 + 7^2 - 7^2 sn^ + k^ 

This is the Lame equation for a particular value of the coefficients that make it solvable in 
terms of Jacobi functions |^. We summarize here the results for the mode functions. The 
derivations are given in ref. UlSj . 

Since the coefficients of eq.( [4.2|) are periodic with period 2u!, the mode functions can be 
chosen to be quasi-periodic (Floquet type) with quasi-period 2u. 

17,(r + 2o^) = e*^('')t/,(r), (4.3) 

where the Floquet indices F{q) are independent of r. In the allowed zones, F{q) is a 
real number and the functions are bounded with a constant maximum amplitude. In the 
forbidden zones F{q) has a non-zero imaginary part and the amplitude of the solutions either 
grows or decreases exponentially. 

Obviously, the Floquet modes Dg(r) cannot obey in general the initial conditions given 
by (^.3|) and the proper mode functions with these initial conditions will be obtained as 
linear combinations of the Floquet solutions. We normalize the Floquet solutions as 


‘P,(t) = 0 . (4.2) 


(7,(0) = 1 . 


We choose Uqir) and Uq{—T) as an independent set of solutions of the second order differ¬ 
ential equation ( [4.2|) . It follows from eq. (|4.3|) that Uq{—T) has —F{q) as its Floquet index. 

We can now express the modes (Pq^r) with the proper boundary conditions [see eq.( |3.3| )] 
as the following linear combinations of fJq(r) and Uq{—T) 


Vqij) 



Uqi-r) + 



Uq{T) 


(4.4) 
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where Wg is the Wronskian of the two Floquet solutions 


W, = ^ -2U,(0) . 


Eq . (14.21) corresponds to a Schrodinger-like equation with a one-zone potential |^. We 


hnd two allowed bands and two forbidden bands. The allowed bands correspond to 


-1 - ^ < < 0 and ^ < +00 , 

and the forbidden bands to 

-oo < < -1 - ^ and 0 < < ^ . 


The last forbidden band is for positive and hence will contribute to the exponential growth 
of the fluctuation function S(r). 

The mode functions can be written explicitly in terms of Jacobi "d-functions for each 
band. We hnd for the forbidden band, 


jj — g-T a/i+»?o z{2K{k)v) ^ 4 ( 0 ) 'di(n + 2 ui) 


(4.5) 


where n is a function of g in the forbidden band 0 < g < ^ dehned by 

« = cii(2K(k) v,k) , 0 < II < i 


(4.6) 


and Z{u) is the Jacobi zeta function [^. It can be expanded in series as follows 

00 

2 K{k) Z{2K{k) v) = Att - — sin(2?7,7rn) (4.7) 

n=i 1 ~ 7 "■ 

where g = The Jacobi 'd-functions can be expanded in series as follows 


§^{v\q) = 2 ^ (-1)^^+^ qin-i/2f gi^^2n - l)7rn , 

n=l 

00 

■di{v\q) = 1 -I- 2 ^ ( — 1)"' cos(2n7rn) . 


n=l 


We explicitly see in eq.( [4.5| ) that Uq{T) factorizes into a real exponential with an exponent 
linear in r and an antiperiodic function of r with period 2a;. Recall that 

'di{x + 1) =, 'd4,{x + 1) =+'d4{x) . (4.8) 

We see that the solution Ugir) decreases with r. The other independent solution Uq{—T) 
grows with r. 

The Floquet indices can be read comparing eq. (^^) , (|4.5|) and (|4.8|) , 

F(g) = 2i K{k) Z{2K{k) v)±tt . 
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Uq{T) turns out to be a real function in the forbidden band. It has real zeroes at 

T = 2uj{n — v) , n eZ . 

and complex poles at 

r = 2u;ni + (2^2 + l)a;' , ni, n 2 eZ . (4.9) 

where oj' is the complex period of the Jacobi functions. Notice that the pole positions 
are g-independent, and that Uq{T) becomes an antiperiodic function on the borders of this 


Uq{T)\q=o = cn(r-\/lT^, k) 

[vUqir)] = sn(r^lT^, k) , 


forbidden band, g = 0 and g = %. We hnd using eq.([4.5[) and ref. 129 


respectively. 

The functions Uq{T) transform under complex conjugation in the forbidden band as 


[I'.M]' = u,(t) . 


(4.10) 


For the allowed band ^ < q < oo, we hnd for the mode functions 


U,(t) 


= e 


2lj n K{k) 


MO) v + 




(4.11) 


where 


q=^^iZ — (2K'(k)v,k') , 
* sn 


o <.<2 


^ Vo 
oo>,>^ 


We see that Uq{T) in this allowed band factorizes into a phase proportional to r and a 
complex periodic function with period 2a;. This function 17g(r) has no real zeroes in r 
except when g is at the lower border q = ^ ■ Its poles in r are g-independent and they are 
the same as those in the forbidden band [see eq. (^^) ]. 

The Floquet indices can be read off by comparing eq. (|4.3|) , (|4.8|) and ([4.11|) 


These indices are real in the allowed band. 

The functions Uq{T) transform under complex conjugation in the allowed band as 


|fi,a)i* = u ,(- t ) . 


(4.12) 
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Obviously these modes will give contributions to the fluctuation S(r) which are always 
bounded in time and at long times will be subdominant with respect to the contributions 
of the modes in the forbidden band that grow exponentially. 

The form of these functions is rather complicated, and it is useful to And convenient 
approximations of them for calculational convenience. 

The expansion of the "d-functions in powers of g = e~'^^ ik)/K{k) converges quite rapidly 
in our case. Since 0 < k < 1/a/2 [see eq.(|4.1|)], we have 


0 < g < e“^ = 0.0432139.... 
q can be computed with high precision from the series [^ 

g = A + 2 A® + 15 + 150 A^=^ + 1707 A^^ + ... , 


where (not to be confused with the coupling constant) 


1 1 - \/F 

2 l + y/W' 


We And from 


eq. (|0|) 


2 (1 + riy/* + (1 + rtoV2)V‘ ' 


The quantity A can be computed and is a small number: for 0 < 770 ^ 00 , we And 0 < A < 
0.0432136.... Therefore, to very good approximation, with an error smaller than ~ 10“^, 
we may use: 


. 1 (1 + riiy/^ - (1 + ril/iyf* 

’ 2(1+ ,g)i/4 + (1 + ^l/2y/i ■ 

We And in the forbidden band from eq. (|4.5|) and 


(4.13) 


JJ — g-4T q sin(27ra) [l+2q (cos 27TO-2)+0 (i}2)] 


l-2g sin7r(n + ^) 


1 — 2g cos(^) sinTTO 


1 + 0(g' 


(4.14) 


where now we can relate n to g in the simpler form 

1 — 4g sin^ ttu + 4g^ sin^ ttu (1 + 4 cos^ vrn) + O(g^) 


ho 

g = -j= cos Tiv 
V2 


which makes it more convenient to write q{v) in the integrals, and 


^ = \/l + [1 - 4 g + 125^ + 0(9*) 


( 4 . 16 ) 


where 0 < n < 

The Floquet indices can now be written in a very compact form amenable for analytical 
estimates 


F(g) = Ain q sin(27rn) 


1 + 2 g cos 2nv + O(g^) 


+ vr . 
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In this approximation the zero mode (^4.1|) becomes 


ti{t) = rjo cos 1 - 4g sin2(^) + 0{q‘^ 


(4.16) 


This expression is very illuminating, because we hnd that a Mathieu equation approxi¬ 
mation, based on the hrst term of eq.( t4.16| ) to the evolution of the mode functions is never 
a good approximation. The reason for this is that the second and higher order terms are 
of the same order as the secular terms in the solution which after resummation lead to the 
identihcation of the unstable bands. In fact, whereas the Mathieu equation has infinitely 
many forbidden bands, the exact equation has only one forbidden band. Even for small g, 
the Mathieu equation is not a good approximation to the Lame equation p8 . 


From eq. (|4.11|) analogous formulae can be obtained for the allowed band 


Ufir) = e 


coth 

2u 


TrK'(k) 

K(k) 


1 — 2q 1 — 2q cos ^ — 2iv log q 


1 — 2q cos(—) 1 — 2g cosh(2n logg) 


l + Oif) 


where 


q = 


VVo 


'ho^ + 2' 


1/4 


23/2 


■K’ik) 

K(k) 


|l -I- 2g cosh(2r> logg) -|- 0(g^)| 


Here, 


0 < n < 


oo > g > 


Vo 

V2- 


Note that eq. (|4.15|) holds in all bands. 

We can now estimate the size and growth of the quantum fluctuations, at least for 
relatively short times and weak couplings. For small times 0 < r < Ti (to be determined 
consistently later) and small coupling g « 1, we can safely neglect the back-reaction term 
gS(r) in eq. (|3.2|) and express the modes (Pg^r) in terms of the functions Ugi^r) and Uq{—r) 
for this however, we need the Wronskian, which in the forbidden band is found to be given 
by: 

1 d ^ 'di(n) „ r. -;t cn dn 

uj dv 'dfiv) ’ sn 


W, = 


log = -2^1 + ril ^^^{2vK{k), k) . 


In terms of the variable g^ this becomes, after using eq. (|4.6|) : 


Wg 


-2q 




f + 1 + g^ 

2a _ 

2 y 


This Wronskian is regular and non-zero except at the four borders of the bands. 

We hnd from eq. (|4.4|) that |(pg(r)p is given by 

1 ( AtD? I 

f,ir) ||f/,(r) + U,(-r)f + ^ (4.17) 
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where we took into account eqs. (|4.10| and ( |4.12|) . Notice that both terms in the rhs of 


eq.( 4.17 ) are real and positive for real q. For very weak coupling and after renormalization, 
the contribution to 5 'S(r) from the stable bands will always be perturbatively small, while the 
contribution from the modes in the unstable band will grow in time exponentially, eventually 
yielding a non-perturbatively large contribution. Thus these are the only important modes 
for the fluctuations and the back-reaction. An estimate of the preheating time scale can be 
obtained by looking for the time when 5 'Fl(r) is of the same order of the classical contributions 
to the equations of motion. In order to obtain an estimate for the latter, we consider the 
average over a period of the classical zero mode: 


1 +< v \ t ) >= (l + p2) 


2E{k) 

K{k) 


- 1 


which yields for small and large initial amplitudes the following results 

< y(r) >’*=•“ I 

0.4569 ... . 

Therefore the average over a period of is to a very good approximation r]Q/2 for all 

initial amplitudes. This result provides an estimate for the preheating time scale ri; this 
occurs when 5 'S(ri) (1 -|- ril/2). Furthermore, at long times (but before (1 -|- rjQ/2)) 

we need only keep the exponentially growing modes and (y'S(r) can be approximated by 


gT.est{r) = ^ q^dq 
4 JO 


1 


<? L 


1 H- - 

g 


Un 


-T 


Moreover, choosing r such that the oscillatory factors in Uq{—T) attain the value 1 (the 


envelope), and using eg. ([4-51) we hnally obtain: 

1 - '^n 

i r 

q dq- 


Ea cV- — 


est—env 


1 


1 

a 


Akll 

1 H- - 

W2 

''q 


qiT a/T+TI Z{2K{k)v,k) 


(4.18) 


where v depends on the integration variable through eq.(|4.6|) 


The Jacobi Z function can be accurately represented using eq. (|4.7|) 

Z{2K{k) v,k) = Aq sin 2nv [1 — 2 g (2 — cos 2?™)] -|- 0{q^) . 
where we recall that q < 0.0433. 

The integral ( |4.18|) will be dominated by the point q that maximizes the coefficient of r 
in the exponent. This happens at g = gi, v = Vi, where 


= 2 ^0 (1 - q) + 0{f) 

Z{2K{k) ui, k) = 4g (1 - 4g) + 0{f) 

We can compute the integral (|4.18|) by saddle point approximation to hnd: 


(4.19) 
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est—env 


,r = 



r. 4 ^ 2 1 


1 1 n 

. 


2 n 


^ST^yT^q (l-4(f) 


91 


f dq ^ ^ 0 ' (i-6« [1 + 0 (g)] 

J —CO 

(1 + g) 




l + !% 


W2 


38 T -y/l+ryg q (l-4q) 


l + 0(-) 
r 


64 (l + r/g)V4 

We can relate q to rjo using eq. ([4.13|) , and we have used the small q expansion 
d^Z 


dv"^ 


{2K{k) ui, k) = —IGvr g (1 — 4g) + 0(g'^) 


dq I 


hoTT 


(1 + g) + 0(g2) . 


In summary, during the preheating time where parametric resonance is important, 
T,est-env{T) cau be represented to a very good approximation by the formula 


Zipaf — 


est—env 


,r = 






(4.20) 


where B and N are functions of t]o given by 

B = 8 \/l + Vo 9 (1 - 4g) + 0{f) , 
64 (l + r^2)i/4y^^^ 


N = 



r, 402 ] 

v-^ 

ho 

1 1 11 

^ W2 

. 



^1/2 

f 

4 /7 (4 + 3 7]^) 


TT 


ho (1 + 


[l + 0(g)] 


(4.21) 


and eg. ([4.131) gives g as a function of r]Q. This is one of the main results of ref. 


We display in Table I below some relevant values of g, B and N as functions of t]q. 

We notice that the limiting values of B and N for r]Q ^ 00 yield a very good approxi¬ 
mation even for ~ 1- Namely, 


E(t) 

with the asymptotic values given by 


T N'oo 


(4.22) 


Boo = 8e-"(l - 4e-")[l + O(77o“ )] = 0.285953 ... [1 + 0{r]Q^)] 


12 


Noo = + Oivf)] = 3.147... [1 + Oi'qf)]. 

\/'ir 


(4.23) 


These rather simple expressions ( [4.20| - ^?23D allow us to perform analytic estimates with 
great accuracy and constitute one of our main analytic results. The accuracy of this result 
will be discussed below in connection with the full numerical analysis including back-reaction. 
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Using this estimate for the back-reaction term, we can now estimate the value of the 
preheating time scale Ti at which the back-reaction becomes comparable to the classical 
terms in the differential equations. Such a time is defined by 5 'S(ri) ~ (1 -|-Pq/ 2). From the 
results presented above, we find 


Ti 


;^log 


N(1 + 1^/2) 
gVB 


(4.24) 


The time interval from r = 0 to r ~ ri is when most of the particle production takes 
place. After r ~ ri the quantum fluctuation become large enough to shut-off the growth of 
the modes and particle production essentially stops. We will compare these results to our 
numerical analysis below. 

We can now use our analytic results to study the different contributions to the energy and 
pressure coming from the zero mode and the quantum fluctuations and begin by analyzing 
the contribution to the energy eo and pressure po from the zero mode pir). 

The dimensionless energy and pressure, (normalized by the factor 2M^/A/j) are given by 
the following expressions. 


«o(r) = i 

P«(r) = 5 


f]'^ + pirf + -^virY 


fY -p{tY - -p{tY 


When the back-reaction term 5 'S(r) can be neglected, we can use eq. (|4.1|) as a good 
approximation to p{t). In this approximation 


eo = -Vo 


1 + 2^0^ 


Po{t) + eo = r^o (l + Vo) snMn^ (^T^Jl + Vo, • 


(4.25) 


The zero mode energy is conserved and the pressure oscillates between plus and minus eo 
with period 2ci;. 

Averaging Po{t) over one period yields 




<po>= — 

zoo Jo 

Inserting eq. (|4.25|) into eq. (|4.26|) yields 


dr po{t) . 


< Po >= --Vo 


1-2^0 +9(l + ^o) 


1 - 


m 

K{k) 


where k is given by eq.(TT). 

< Po > vanishes for small po faster than eo. 


<Po >^°=^ ^Vo + OiVo), 


(4.26) 


(4.27) 
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so that the zero mode contribution to the equation of state is that of dust for small rjo. For 
large t]q we hnd from eq. ( 4.27 ), 


r]o—>oo 4 I 2 

<Po> = Y^Vo + Vo 


2 E{l/V2) 

3 K{l/V2) 


+ 0 ( 1 ) 


where | — | = 0.01437.... The equation of state approaches that of radiation for 

Vo oo: 


rjo->oo 

<Pq> = €o 


0.6092 

ho 


+ 0{Vi 


-4\ 


Thus we see that for small amplitudes the zero mode stress-energy, averaged over an 
oscillation period, behaves as dust while for large amplitudes, the behavior is that of a 
radiation fluid. The ratio < po > /so for zero mode vs.^o is shown in hgure 1. 

The contribution from the /c 7 ^ 0 modes originates in the quantum fluctuations during 
the the stage of parametric amplihcation. 

Since we have fluid behaviour, we can dehne an effective (time-dependent) polytropic 
index 7 (r) as 

7(r).hh + i, 

where renormalized quantities are understood throughout. Within a cosmological setting 

whenever 7 (r) reaches a constant value such equation of state implies a scale factor R{t) = 

2_ 

Rq r 37. 

In the case being studied here, that of Minkowski space, £ is time-independent and hence 
equal to the initial energy density (divided by N and restoring pre-factors) which after a 
suitable choice of the constant C is given by: 


fl , 

“ Ah l2’'» 



As argued before, for weak coupling the important contribution to the quantum fluc¬ 
tuations come from the modes in unstable bands, since these grow exponentially in time 
and give rise to a non-perturbatively large contribution. Thus we concentrate only on these 
modes in calculating the pressure. 

The contribution of the forbidden band to the renormalized p(r) - 1 - £ can be written as 


[p(r) -t- e]unst = 





Ag(^) P 


A<?( 


T] 



After renormalization, the terms that we have neglected in this approximation are pertur- 
batively small (of order g) whereas the terms inside the bracket eventually become of order 
1 (comparable to the tree -level contribution). We now only keep the exponentially growing 
pieces in the mode functions Pqir) and Pgir) since these will dominate the contribution to 
the pressure. This is simplihed considerably by writing to leading order in q 


Ag(^) 
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Averaging over a period of oscillation yields 


[ p ( t ) -|“ ^^unst 


q^dq 1 

A_r [4 jo (47r)^ 2 f2q sin^ TTV 


Ani 

1 H- - 


(4.28) 


^2t i/l+gg Z{2K{k)v,k) 


1 + ^0 + 3 


This integral is similar to the one in eg . (14.181) and we find that they are proportional in 
the saddle point approximation. In fact, 


\p{r) + e]unst = 


2\Mn\ 


4 r 


A 


R 


9^est-env{'l~) ( ^ 


where Eest-env{r) is given by eq. ([4.20|) . 
The effective polytropic index 7 (r) is: 


d(”^) d'^est—envij'') 


12 + 1372 


3^0 Oo+ 2) ■ 

When g'Lest-envi.'Ti) ~ 1 + 7o/2; i-e. at the end of the preheating phase, 7 (r) is given by 

12 + 1372 


7e// OC 


670 


We note here that for very large 70 the effective polytropic index is 7 e// — 13/6 ~ 0(1). It 
is clear then that the physics can be interpreted in terms of two fluids, one the contribution 
from the zero mode and the other from the fluctuations, each with an equation of state that 
is neither that of dust nor of radiation, but described in terms of an effective polytropic 
index. 

We can now use our approximations to obtain an estimate for the number of particles 
produced during the preheating stage. In terms of dimensionless quantities, the particle 
number, defined with respect to the initial Fock vacuum state is given by eq.( p.ll|) . 

This particle number will only obtain a significant contribution from the unstable modes 
in the forbidden band where to leading order in q we can approximate +g(r) and +g(r) by 
its exponentially growing pieces [see eq.( |4.17| )], as follows: 


i+gwr - 


4 a 


4 a 

1 H- - 

>V2 

<? 




\'Pq{T)?' ~ (1 + 


TT n — 


2uj' 


Wq{'T)? + 0{q) 


The total number of produced particles N{t) per volume is given by: 


M = 


N{t) _ r d?q 


Nq{T) . 


The asymptotic behaviour (|3.16|) ensures that this integral converges. 
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Following the same steps as in eq. (^4.18|) and ([4.28|) , we find 


n(t: 


unst 




"^est—env ) 


1 +ho 


+ Vt 


91 


91 


1 ^ + hl 

^A + bril 


idlest —env (r)) . 


where we nsed eq.( |4.19|) and Tjest-enviR) is given by the simple formula ( |4.2CI|) . Notice that 
by the end of the preheating stage, when 5 'Fl(r) ~ 1 + ? 7 q /2 the total number of particles 
produced is non-perturbatively large, both in the amplitude as well as in the coupling 


Mt 


1 (4+|r;2)(l + r;oV2) 


tot 


A 


R 


4 + 577^ 


(4.29) 


The total number of adiabatic particles can also be computed in a similar manner with 
a very similar result insofar as the non-perturbative form in terms of coupling and initial 
amplitude. 


B. Analytic Results in the Hartree and resummed one-loop approximations 


We give in this section the analytic treatment of the Hartree and one-loop equations 
during preheating. The full Hartree equations ( p.l0| )- (p.l4|) in dimensionless variables take 
the form. 


ii + r] + r]^ + 3g E(r) = 0 
^ + 1 + 3 r/(r)2 + 3g E(r) 

</59(0) = 


09(0) = 


h(0) = ho , h(0) = 0 


(4.30) 

^,(r) = 0 , (4.31) 

(4.32) 

(4.33) 

where g H(r) is given by eq.( ^(6D . 

Eqs. (|4.30|) - ([4.33|) only differ from the large N eqs. ( p.l|) - (p.4|) on factors of 3 in some 
coefficients. 

As in sec. IVA, in the weak coupling regime and for times small enough so that the 
quantum fluctuations are not large compared to the ‘tree level’ quantities, we can neglect 
gS(r). Since S(0) = 0, such approximation is expected to hold for an interval 0 < r < ti, 
where ti will be called the preheating time. 

During this interval of time we can then approximate h('r) by the classical solution ( [4.11 ). 
Inserting this elliptic function for ^{t) in eq. (|4.31|) and neglecting hH(r) yields 


dr^ 


+ + 1 + 3 ho cn^ (J^Jl + r|l 


(Pq{T) = 0 


(4.34) 


This is the Lame equation for a particular value of the coefficients that make it solvable in 
terms of Jacobi functions |^. We summarize here the results for the mode functions. The 
derivations are analogous to those given in ref. [0. As for the large N limit, we can choose 
the mode functions here to be quasi-periodic on r (Floquet type). 
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Eq. (|4.34|) corresponds to a Schrodinger-like equation with a two-zone potential . We 
find three allowed bands and three forbidden bands. The allowed bands correspond to 


V3r/o' + 6r/g + 4 + 1 < 


0<g2<^ + 3 


and 

Y^3r/^T”67i^”+4 + l < < -Fco , 

and the forbidden bands to 

-oo <q^< + Qr]l + A + l , < 0 

and 

^ + 3 < g^ < ^377^ + 677^ + 4 + 1 . 

The last forbidden band is for positive q^ and hence will contribute to the exponential growth 
of the fluctuation function S(r). 

The mode functions can be written explicitly in terms of Jacobi "d-functions for each 
band. We find, 


where, 

cLt 

' ■ M±) ' 


(4.35) 


where 0 < n < | is a function of g in the forbidden band ^ + 3 < g^ < + 1 

defined by 


2(ho^ + 1) 


and 


(3{v) = 


sn^ (2K{k) n, k) 

+ 1) 


Vo 


2g2 (g2 - 3)^ 


3(772 + l)"-g2(g2-2) 


cn dn 


2(770 + 1) + ( 1 ?^ — 1 — 770 ) sn^ {2K{k) V, k) sn 


{2K{k)v,k) 


1 7?; 

2K{k)¥i 


(4.36) 


v) ■ 


Eq. (|4.36|) is a third order equation in g^ defining g^ as a function of v. We can express 
its solution in compact form as follows 


where 


and 


g^ = 1 — w{v) [1 + 2 cos (a + 27r/3)] , 
777(7;) = 3(772 + 1) 


1 ^ cn2 (2K{k) 77, k) 


1 

+ 2 


cos 3a = 1 


2 sn^ (2it'(fc) 77, fc) 

9(772 + 1)2 + 3 9(77o2 + 1)2-1 


2777(77)2 


2777 ( 77)3 
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with 0<a<7r/3. At the upper border of the band, = \j3r]Q + br^g + 4 + 1, we have 

V = a = 0. The lower border of the band, ^ + 3, corresponds to n = |, a = tt/3. 
The Floquet index dehned by eq. ( |4.3|) takes here the form 

F{q) = -2iK{k)l3{v) . 

The mode functions in the other bands follow from eq. ([4.35D by analytic continuation in 
V. In addition, explicit and accurate expressions for the Floquet indices as well as the for 
the mode functions can be obtained by expanding in powers of q as in sec. IV A. 

In this approximation valid for times early than the preheating time, the Hartree and 
the resummed one-loop approximation are indeed identical. 


C. Numerical Results 


We now evolve our equations for the zero and non-zero modes numerically, including the 
effects of back-reaction. We will see that up to the preheating time, our analytic results 
agree extremely well with the full numerical evolution. 

The procedure used was to solve equations (|3.1|, |^) with the initial conditions (|373| , 
3^ , |3.5| ) and (|3.6| ) using a fourth order Runge-Kutta algorithm for the differential equation 
and an 11-point Newton-Cotes integrator to compute the fluctuation integrals. We tested 
the cutoff sensitivity by running our code for cutoffs A/|M/j| = 100, 70, 50, 20 and for very 
small couplings (which is the case of interest. We found no appreciable cutoff dependence. 
The typical numerical error both in the differential equations and the integrals are less than 
one part in 10®. 

Figure 2.a shows r 7 (r) vs.r for r^o = 4.0 ,g = 10“^^. For this weak coupling, the effect 
of back-reaction is negligible for a long time, allowing several undamped oscillations of the 
zero mode. Figure 2.b shows (y'S(r) vs. r. It can be seen that the back-reaction becomes 
important when (y'S(r) 1 + rj^/2 as the evolution of '^{t) begins to damp out. This 

happens for r ~ 25 in excellent agreement with the analytic prediction given by eq. (|4.24|) 
Ti = 26.2 ..., the difference between the analytic estimate for Fl(r) given by eq. ([4.20|) and the 
numerical result is less than 1% in the range 0 < r < 30. Figure 2.c shows gA/'(r) vs. r and 
we see that the analytic expression ( |4.29|) gives an approximate estimation AiiMoi 74.6 ... 
for the hnal number of produced particles. 

Figures 2.d-2.f, show gNg^r) for r = 40,120, 200, we see that the prediction of the width 
of the unstable band 0 < g < r/o/v^ is excellent and is valid even for very long times beyond 
the regime of validity of the small time, weak coupling approximation. However, we see that 
the peak becomes higher, narrower and moves towards g ~ 0.5 as time evolves beyond ri. 
This feature persists in all numerical studies of the unbroken phase that we have carried 
out, this changes in the peak width, height and position are clearly a result of back-reaction 
effects. We have searched for unstable bands for 0 < g < 20 and we only found one band 
precisely in the region predicted by the analytic estimate. All throughout the evolution there 
is only one unstable band. The band develops some structure with the height, position and 
width of the peak varying at long times but no other unstable bands develop and the width 
of the band remains constant. For values of g outside the unstable band we find typically 
gNq < 10“^^ at all times. This is a remarkable and unexpected feature. 
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Obviously this is very different from the band structure of a Mathieu equation. The 
Mathieu equation gives rise to an infinite number of narrowing bands, so that quanti¬ 
tative estimates of particle production, etc. using the Mathieu equation approximation 
would be gross misrepresentations of the actual dynamics, with discrepancies that are non- 
perturbatively large when the back-reaction becomes important |^. Since particle produc¬ 
tion essentially happens in the forbidden bands, the quantitative predictions obtained from 
a single forbbiden band and an infinite number, as predicted by WKB or Mathieu equation 
analysis, will yield different physics. 

We have carried the numerical evolution including only the wave-vectors in the unstable 
region and we find that this region of g-wavevectors is the most relevant for the numer¬ 
ics. Even using a cutoff as low as gc = 4 in this case gives results that are numerically 
indistinguishable from those obtained with much larger cutoffs Qc = 70 — 100. The occu¬ 
pation number of modes outside the unstable bands very quickly becomes negligible small 
and for g 4 it is already of the same order of magnitude as the numerical error < 10“^°. 
Clearly this is a feature of the weak coupling case under consideration. Keeping only the 
contribution of the modes in the unstable band, the energy and pressure can be written as 


£ = 


P = 


Qc = 




\rf ^ rf ^ P 
Ah 

0{9)] 

2\Mr\^ 


v\r) -p^ + ^S(r) 


A 

Vo 

V2 


R 




+ l0Wr +v‘^ + Oig) \ -e 


+ 

(4.37) 

(4.38) 


where we have made explicit that we have neglected terms of order g in eqs.( [4.3'iT - |4.38|) . The 
terms multiplied by g in eqs. (|4.37| , [4.38|) become of order 1 during the preheating stage. For 
the parameters used in figures 2, we have checked numerically that the energy ( [f.37|) is con¬ 
served to order g within our numerical error. Figure 2.g shows the pressure 2\Mr\^ p{'t)/^r 
versus r. Initially p(0) = —e (vacuum dominated) but at the end of preheating the equation 
of state becomes almost that of radiation = s/O. 

For very small coupling (g ~ 10“^^), the backreaction shuts-off suddenly the particle 
production at the end of the preheating (see fig. 2c). Later on, (r larger than 100 for 
g ~ 10“^^) the time evolution is periodic in a very good approximation. That is, this non¬ 
linear system exhibits a limiting cycle behaviour. The modulus of the fc-modes do not grow 
in time and no particle production takes place. This tells us that no forbidden bands are 
present for g^ > 0 in the late time regime. 

We have numerically studied several different values of r/o 5 S' finding the same qualitative 
behavior for the evolution of the zero mode, particle production and pressure. In all cases 
we have found remarkable agreement (at most 5% difference) with the analytical predictions 
in the time regime for which 0 < gF(r) < 1. The asymptotic value of the pressure, however, 
only becomes consistent with a radiation dominated case for large initial amplitudes. For 
smaller amplitudes po = 1 we find that asymptotically the polytropic index is smaller than 
4/3. This asymptotic behavior is beyond the regime of validity of the approximations in 
the analytic treatment and must be studied numerically. This polytropic index depends 
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crucially on the band structure because most of the contribution comes from the unstable 
modes. 


In ref. ITHI results of ref. IHI are rederived with a different renormalization scheme. 


V. THE BROKEN SYMMETRY CASE 
A. Analytic Results 

As in the unbroken case, for << 1 we can neglect 5 'S(r) in eq. (p.7|) until a time T 2 
at which point the fluctuations have grown to be comparable to the ‘tree level’ terms. The 
zero mode equation then becomes, 


7] — T] + T]^ = 0 . 

which correspond to the evolution on the classical potential 

^ = (5.1) 

with the initial conditions ( |3.4|) . We then hnd for 0 < ?7o < 1, 



k 


1-ho 



(5.2) 


Notice that //(r) 

For 1 < rjQ < 

? 7 (r) = r]o dn (Tr]o/V2, fc) 

k = ^2{l-7jo-^) . (5.3) 

This solution follows by shifting eq.( b.2|) by a half-period and changing rj^ ^ 2 — 7]^. It has 
a period 2a; = ^^K(k). For r]Q —> 1, 2a; —>■ TTy/2 and the oscillation amplitude vanishes, 
since = 1 is a minimum of the classical potential. 

For po > we obtain 


has period 2a; = 
^J2 we hnd 


i_2a 
^ 2 


The elliptic modulus k is given by eq.(|5.2|). 


r7(r) 


k 


ho cn 



ho 

pw - 1 ) 


(5.4) 


This solution has 4a; = as period. 

V^o-i 
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The solutions for rjQ < a /2 and rjQ > a /2 are qualitatively different since in the second case 
//(r) oscillates over the two minima rj = ±1. In the limiting case r/o = these solutions 
degenerate into the instanton solution 


V{r) = 


V2 

coshr ’ 


and the period becomes inhnite. 


Inserting this form for /^(r) in eq. (|3.8|) and neglecting 5 'S(r) yields for 0 < r/o < 1, 




1 + 


do 




ip,{T) = 0 . 


( 6 . 6 ) 


This is again a Lame equation for a one-zone potential and can also be solved in closed form 
in terms of Jacobi functions. We summarize here the results for the mode functions, with 


the derivations again given in ref. [15 


As for unbroken symmetry case, there are two allowed bands and two forbidden bands 
The allowed bands for 0 < < 1 correspond to 


0 < < — and 1 — — < < J-cx) , 

_y_2 2 - 


do 


and the forbidden bands to 


oo<g2<0 and < g^ < 1 _ . 

-H - 2 -'^- 2 


(5.6) 


The last forbidden band exists for positive g^ and hence contributes to the growth of S(t). 
The Floquet solutions obey eqs. ( [4.3|) and the modes ^qir) can be expressed in terms of 


t/g(r) and Uq{—r) following eq. (|4.4|) . 

It is useful to write the solution I 

2 2 _ 

band ^ ^ ~ ^ after some calculation (see ref. |T5[]), 


It is useful to write the solution f/q(r) in terms of Jacobi "d-functions. For the forbidden 
2 ^ ^2 


Uqir) = e 


\J 1 -^ Z{2K{k)v) '^ 3 ( 0 ) '^2{V + 2 /) 


(5.7) 


where 0 < u < ^ is related with g through 


= 1 - Y - (1-do) (2 77(fc)u,fc) , 


and fc is a function of rjo as defined by eq. (|5.2|) . 

We see explicitly here that Uq{T) factorizes into a real exponential with an exponent 
linear in r and an antiperiodic function of r with period 2a;. 

The Floquet indices for this forbidden band are given by 


F(g) = 2 i K{k) Z{2K{k) u) ± tt . 
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For the allowed band 1 — % < 5'^ < +cxo , we find for the modes, 


T 1 \ 

U„{t) = e 


^ 3 ( 0 ) 


2 a; 


urn Ui) ’ 


where q and a are related by: 


q. = 


1 _ !la 
2 




sn a 


with > a > 0. The Floquet indices for this first allowed band are given by 

Vi-^o 

™>-s© 

2 

Analogous expressions hold in the other allowed band, 0 < < y- 

U M 1(0 ’^3(0) U'^) 

^ MB Mi) ■ 

Here, g = ^ sn {a\Jl — k'^ and > 0, and the Floquet indices for this band 

are given by 

For r^o ~ 1 the situation is very similar to the unbroken symmetry case; the zero mode 
oscillates quasi-periodically around the minimum of the tree level potential. There are effects 
from the curvature of the potential, but the dynamics can be analyzed in the same manner 
as in the unbroken case, with similar conclusions and will not be repeated here. 

The case rjQ << 1 is especially interesting P,p!5| for broken symmetry because of new and 
interesting phenomena that has been recently associated with symmetry restoration 

In this limit, the elliptic modulus k [see eq. ( b.2|) ] approaches unity and the (real) period 
2a; grows as 

2a; ~ 2K{k) + O('go) — 2 In + ^(^o) 

In this limit, both the potential in eq. (|5.5|) and the mode functions ( b.7|) can be approximated 
by hyperbolic functions ||29|| : 


dn(r-^l - f, k) 


= coshr + 0{r]l) 


u 


Z{u) = tanhw “ ^ + ^(^o) 


where 
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coshw = 


A = log 



(5.8) 


Using the imaginary Jacobi transformation p9|, 




Kik) _zLi£W^2 .Kik) ... 


■KK{k ^) 


rrK{k) 


where q = e , q = e and the series expansions ^ 


■(93(^1^) = 1 + 2 ^ cos(27r?w) 

n=\ 

di{v\q) = + \\q) ^ 

we can derive expressions for the mode functions Uj(r) valid for small rjQ-. 


Uq{r) = e 


tanhii^l+^^ I — ^ coshw COsh(M + 2t) 


1 —^ coshti 


l+OhS 


(5.9) 


Here u is related with q through eq. (|5.8|) . We see that the function Ug{—T) grows with r 
almost as e’’ for q near the lower border of the forbidden band m ~ A. This fast growth 
can be interpreted as the joint effect of the non-periodic exponential factor in eq.( b-7|) and 
the growth of the periodic ■d-functions. Since the real period is here of the order A, the 
two effects cannot be separated. The unstable growth for t < uj also reflects the spinodal 
instabilities associated with phase separation |^. 

In this case, there is a range of parameters for which the quantum fluctuations grow 
to become comparable to the tree level contribution within just one or very few periods. 
The expression (|5.9|) determines that S(r) from the contributions of modes near the 

lower edge of the band. The condition for the quantum fluctuations to become of order 
1 within just one period of the elliptic function is ~ 1 which leads to the conclusion 
that for 77 ( 0 ) < the quantum fluctuations grow very large before the zero mode can 
actually execute a single oscillation. In such a situation an analysis in terms of Floquet 
(quasi-periodic) solutions is not correct because the back reaction prevents the zero mode 
from oscillating enough times for periodicity to be a reasonable approximation. 

We now analyze the behavior of the pressure for the zero mode to compare to the previous 
case. In the approximation where eqs. (^.2| ), (|5.3|) and ( |5.4| ) hold and adjusting the constant 
C in the dehnition of the energy, we have 


to = J {nl - 1)" . 
Po(’') = -to + 


(5.10) 


Inserting eqs. (|5.2|) , (|5.3|) and (|5.4| ) in eq.( b.lOI) yields 
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(5.11) 


0 < ?7o < 1 : Po(r) = -eo 


1 — 8 sn^cn^ [{t + K)\J1 — 


1 < < V2 : Pq{t) = -eo 1-8 sn^cn^ fc) 


^0 


> ^2 : po('r) = -eo 


1 — 8 sn^dn^ — 1, 


Notice that the functional form of the elliptic modulus fc as a function of r/o is different in 
each interval [see eqs. (|5.2|) , (|5.3|) and ( ^.41) ]. 

Let us now average the pressure over a period as in eq.([4.26|). We find 


< Po > = eo 


-2 


K{k ); p 


-1 


for 0 < po < , 


<Po> =^0^2 


2^ 


(1-2A;") 1- 


K{k)) 


- 1 


for po > . 


(5.12) 


The dimensionless energy eo tends to 1/4 both as r/o ^ 0 and r/o ^ in both cases we 
hnd using eq. (|5.12|) 


< Po > 
eo 


-1 - 


16 


3 log|i -eol 


+ O(eo--). 


This is result is recognized as vacuum behaviour in this limit. 
For po —^ 1, eq. (|5.12|) yields, 


0 (r,„ - 1)2 , 

eo 

That is a dust type behaviour, which is consistent with the small amplitude limit of the 
unbroken symmetry case studied before. 

Finally, for rj^ ^ oo, when the zero mode is released from high up the potential hill, we 
find that the pressure approaches radiation behaviour (from above) 


< PO ^ 170^00 1 ^ 

eo “ 3 3 



4 + 0(4) 

m Vo 


1 

3 


0.86526. 

Vo 


0 ( 4 ) 

Vo 


Figure 3 shows < po > /eq vs. £o- 

As mentioned before, we expect that for po << 1 the conclusion will be modified dra¬ 
matically by the quantum corrections. 
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B. Numerical Results 


In the region ?7o ~ 1 the analytic estimates are a good approximation for large times and 
weak couplings. We have studied numerically many different cases with rjo > 0.5 and weak 
coupling and conhrmed the validity of the analytic estimates. These cases are qualitatively 
similar to the unbroken symmetry case with almost undamped oscillations for a long time 
compatible with the weak coupling approximation and when 5 'S(r) grows by parametric 
amplification to be of order one with a consequently large number of produced particles and 
the evolution of the zero mode damps out. 

However as argued above, for t]q << 1 the analytic approximation will not be very 
reliable because the quantum fluctuations grow on a time scale of a period or so (depending 
on the coupling) and the back-reaction term cannot be neglected. Thus this region needs to 
be studied numerically. 

We numerically solved equations (|3.7|) -( ^78|) with the initial conditions (p.9|) , (|3.10|) and 
( p.lOj ). The numerical routines are the same as in the unbroken symmetry case. Again we 
tested cutoffs A/\Mji\ = 100,70,50,20 and for very small couplings (which is the case of 
interest, g = 10“®-- -g = 10“^^) we found no appreciable cutoff dependence, with results 
that are numerically indistinguishable even for cutoffs as small as gc ~ 2 . The typical 
numerical error both in the differential equations and the integrals are the same as in the 
unbroken case, less than one part in 10®. 

We begin the numerical study by considering first the case of very small coupling and 
gQ << 1; later we will deal with the case of larger couplings and initial values of the zero 
mode. Figure 4.a shows g{T) vs. r for r/o = 10“®, g = 10“^®. In this case we see that 
within one period of the classical evolution of the zero mode, gFl(r) becomes of order one, 
the quantum fluctuations become non-perturbatively large and the approximation valid for 
early times and weak couplings breaks down. Fig. 4.b shows ^^(r) and £g. 4.c shows gM{T) 
vs. r for these parameters. We hnd that only the wave vectors in the region 0 < g < 1 are 
important i.e. there is only one unstable band whose width remains constant in time. This is 
seen in figs. 4.d-f, which show the particle number (defined with respect to the initial state) 
as a function of wave vector for different times, gNq^r = 30) ,gNq{T = 90) ,gNq{T = 150) 
respectively. Although the analytic approximation breaks down, the prediction eq. (^.6|) for 
the band width agrees remarkably well with the numerical result. As in the unbroken case, 
the band develops structure but its width is constant throughout the evolution. As can 
be seen in these figures the peak of the distribution becomes higher, narrower and moves 
towards smaller values of g. The concentration of particles at very low momentum is a 
consequence of the excitations being effectively massless in the broken symmetry case. The 
features are very distinct from the unbroken symmetry case, in which the peak approaches 
g ^ 0.5. 

We found in all cases that the asymptotic behavior corresponds to 


Al^(r) = -1 + ■g^(r) -h gF(r) 


lim -7 


0 


This is a consistent asymptotic solution that describes massless “pions” and broken symme¬ 
try in the case g{oo) ^ 0. 

For times r ~ 100 — 150 the value of the zero mode is somewhat larger than the initial 
value: //(r = 150) k, 2 x 10“®. This result, when combined with the result that the average 
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of the effective mass approaches zero is clearly an indication that the symmetry is broken. 
We fonnd nnmerically that the hnal valne of the zero mode depends on the initial valne and 
the conpling and we will provide nnmerical evidence for this behavior below. 

Figure 4.a, presents a puzzle. Since the zero mode begins very close to the origin with 
zero derivative and ends up very close to the origin with zero derivative, the classical energy 
of the zero mode is conserved. At the same time, however, the dynamical evolution results in 
copious particle production as can be seen from £g. 4.c. We have shown in a previous section 
that the total energy is conserved and this was numerically checked within the numerical 
error. Thus the puzzle: how is it possible to conserve the total energy, conserve the classical 
zero mode energy and at the same time create 0(1/g) particles? The answer is that there is 
a new term in the total energy that acts as a ‘zero point energy’ that diminishes during the 
evolution and thus maintains total energy conservation with particle production. The most 
important contribution to the energy arises from the zero mode and the unstable modes 
0 < q < Qu- The energy and pressure are given by (adjusting the constant C such that the 
energy coincides with the classical value): 


^ ^ {^cl £n -k Sc d(fi')} 

£Ar(r) = 2g q^dq Qg Ag(r) 

^ 0 


£c{t) = -s(r) 


-I - + M\r) - 


1 I 2 


p{t) = 


^ lo 


||^,(r)p + |0(r)p 


+ f] +0{g) 


— e 


(5.13) 

(5.14) 

(5.15) 


M (t) = -1 + r] (t) + 5 'S(r) 


where Al(r)^ is the effective squared mass of the — 1 “pious” and again 0{g) stand for 
perturbatively small terms of order g. The terms displayed in (|5.13| - |5.15|) are all of 0{1) 
during the preheating stage. 

We hnd that whereas Sn^t) grows with time, the term ec('r) becomes negative and 
decreases. In all the cases that we studied, the effective mass A4 (r) approaches zero asymp¬ 
totically; this is seen in figure 4.g for the same values of the parameters as in figs. 4.a-c. This 
behavior and an asymptotic value r/oo 7 ^ 0 is consistent with broken symmetry and massless 


pions by Goldstone’s theorem. The term Sc(t) in eq.( b.l4| ) can be identihed with the ‘zero’ 
of energy. It contributes to the equation of state as a vacuum contribution, that is pc = —Sc 
and becomes negative in the broken symmetry state. It is this term that compensates for 
the contribution to the energy from particle production. 

This situation is generic for the cases of interest for which r^o << 1, such is the case 
for the slow roll scenario in inflationary cosmology. Figs. 4.h-k show Sdir) vs. r, e]d{T) 
vs. r, Scir) vs. r and 2 |M^p P('^) same values of parameters as £g. 4.a. 

The pressure has a remarkable behavior. It begins with p = —e corresponding to vacuum 
domination and ends asymptotically with a radiation-like equation of state p = e/3. A simple 
explanation for radiation-like behavior would be that the equation of state is dominated by 
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the quantum fluctuations which as argued above correspond to massless pions and therefore 
ultrarelativistic. It must me noticed that we obtain a radiation-like equation of state in spite 
of the fact that the the distribution is out of equilibrium and far from thermal as can be 
seen from figs. 4.d-f. 

An important question to address at this point is: why does the zero mode reach an 
asymptotic value different from the minimum of the effective potential? The answer to 
this question is that once there is profuse particle production, the zero mode evolves in a 
non-equilibrium bath of these excitations. 

Through the time evolution, more of these particles are produced and the zero mode 
evolves in a highly excited, out of thermal equilibrium state. Furthermore we have seen 
in detail that this mechanism of particle production modifies dramatically the zero point 
origin of energy through the large and negative term Sc{t) and therefore the minimum of the 
effective action, which is the appropriate concept to use for a time dependent evolution as 
the present one. The final value reached by the zero mode in the evolution will be determined 
by all of these non-perturbative processes, and only a full non-linear study (including the 
back-reaction of modes on themselves) captures the relevant aspects. As we have argued 
above, approximations based on Mathieu-type equations or the WKB approximation are 
bound to miss such important non-linear processes and will lead to an incomplete picture 
of the evolution. 

These time dependent processes cannot be studied using the effective potential. For 
example, the profuse particle production taking place here is a feature completely missed 
by the effective potential. We And that the effective potential is an irrelevant quantity to 
study the dynamics 


VI. BROKEN SYMMETRY AND ITS QUANTUM RESTORATION AT 

PREHEATING 


The numerical result depicted by figures 4.a-c, [|,|^], has motivated the suggestion that 
the growth of quantum fluctuations is so strong that the non-equilibrium fluctuations restore 
the symmetry aa- The argument is that the non-equilibrium fluctuations given by the 
term g F(r) in eq. (p.8|) for the mode functions grow exponentially and eventually this 
term overcomes the term — |m^| leading to an effective potential with a positive mass squared 
for the zero mode [see eq. O). 

Although this is a very interesting suggestion, it is not borne out by our numerical 
investigation for rjQ < 1. The signal for broken or restored symmetry is the final value 
of the zero mode when the system reaches an equilibrium situation. Any argument about 
symmetry restoration based solely on the dynamics of the fluctuation term g t]{t) S(r) is 
incomplete if it does not address the dynamics of the zero mode. In particular for the case 
of figures 4.a-c, the initial value of the zero mode Pq ^ 0 and the final value is very close to 
the initial value but still different from zero. 

At the same time, the asymptotic effective mass of the ‘pions’ is on average zero. Clearly 
this is a signal for symmetry breaking. Because the initial and final values of the order 
parameter are so small on the scale depicted in the figures, one could be tempted to conclude 
that the symmetry originally broken by a very small value of the order parameter is restored 
asymptotically by the growth of non-equilibrium fluctuations. To settle this issue we show a 
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different set of parameters in figures 5.a,b that clearly show that the final value of the order 
parameter r^oo 7 ^ 0, while the effective mass of the pions M(t) — >■ 0. Here tjq = 0.01 and 
g = 10“^, and asymptotically we hnd r]{T = 150) ~ 0.06, the average of the effective mass 
squared = 0 and the symmetry is broken, despite the fact that the fluctuations have 

grown exponentially and a number of particles 0(1/g) has been produced. 

The reason that the symmetry is not restored is that when the effective mass becomes 
positive, the instabilities shut off and the quantum fluctuations become small. When this 
happens gT, is no longer of order one and the instabilities appear again, producing the oscil¬ 
latory behavior that is seen in the hgures for 5 'H(r) at long times, such that the contributions 
of the oscillatory terms average to zero. It is rather straightforward to see that there is a 
self-consistent solution of the equations of motion for the zero mode and the fluctuations 
with constant rjoc and Al^(oo) = 0. Eq. ( |3.7| ) takes the asymptotic form ^ 


hoo 


-^ + via + 9 ^{ 00 ) 


0 . 


In addition, eq. (|^) yields when Al(cxo)^ = 0, 

(p,(r) + B, e*"" , 


where Ag and Bg depend on the initial conditions and g. We get from eqs. (|3.6|) and (|3.11|) , 

n 


= 2 dQ-g S{vo) , 

JO + 


( 6 . 1 ) 


where S{go) = 1(1- gl) [log 




+ 7(1 + Vo) 


ArgTh 


l-^n 


arctg 




We see that the value of g^o depends on the initial conditions. Whereas the last term 
in eq. (|6.1|) is perturbatively small, the contribution from the produced particles is non- 
perturbatively large, as Ng{oo) ^ 1/g for the unstable wavevectors. Thus the asymptotic 
value of the zero mode is drastically modihed from the tree level v.e.v (in terms of renor¬ 
malized parameters) because of the profuse particle production due to the non-equilibrium 
growth of fluctuations. 

Another way to argue that the symmetry is indeed broken in the hnal state is to realize 
that the distribution of “pions” at late times will be different than the distribution of the 
quanta generated by the fluctuations in the a held, if for no other reason than that the pions 
are asymptotically massless while the a quanta are massive, as long as g^o is different from 
zero. If the symmetry were restored during preheating, these distributions would have to be 
identical. 

In the situation of ‘chaotic initial conditions’ but with a broken symmetry tree level 
potential, the issue of symmetry breaking is more subtle. In this case the zero mode is 
initially displaced with a large amplitude and very high in the potential hill. The total 
energy density is non-perturbatively large. Classically the zero mode will undergo oscillatory 
behavior between the two classical turning points, of very large amplitude and the dynamics 
will probe both broken symmetry states. Even at the classical level the symmetry is respected 
by the dynamics in the sense that the time evolution of the zero mode samples equally both 
vacua. This is not the situation that is envisaged in usual symmetry breaking scenarios. For 
broken symmetry situations there are no hnite energy held conhgurations that can sample 
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both vacua. In the case under consideration with the zero mode of the scalar held with very 
large amplitude and with an energy density much larger than the top of the potential hill, 
there is enough energy in the system to sample both vacua. (The energy is proportional to 
the spatial volume). Parametric amplihcation transfers energy from the zero mode to the 
quantum huctuations. Even when only a fraction of the energy of the zero mode is transferred 
thus creating a non-perturbatively large number of particles, the energy in the huctuations 
is very large, and the equal time two-point correlation function is non-perturbatively large 
and the held huctuations are large enough to sample both vacua. The evolution of the zero 
mode is damped because of this transfer of energy, but in most generic situations it does 
not reach an asymptotic time-independent value, but oscillates around zero, sampling the 
tree level minima with equal probability. This situation is reminiscent of hnite temperature 
in which case the energy density is hnite and above a critical temperature the ensemble 
averages sample both tree level vacua with equal probability thus restoring the symmetry. 
In the dynamical case, the “symmetry restoration” is just a consequence of the fact that 
there is a very large energy density in the initial state, much larger than the top of the tree 
level potential, thus under the dynamical evolution the system samples both vacua equally. 
This statement is simply the dynamical equivalent of the equilibrium hnite temperature 
statement that the energy in the quantum huctuations is large enough that the huctuations 
can actually sample both vacua with equal probability. 

Thus the criterion for symmetry restoration when the tree level potential allows for 
broken symmetry states is that the energy density in the initial state be larger than the 
top of the tree level potential. That is when the amplitude of the zero mode is such that 
I/(po) > 1^(0). In this case the dynamics will be very similar to the unbroken symmetry 
case, the amplitude of the zero mode will damp out, transferring energy to the quantum 
huctuations via parametric amplihcation, but asymptotically oscillating around zero with a 
fairly large amplitude. 

To illustrate this point clearly, we plot in hgs. 6 and 7, r]{t) and E(f) for r^o = 1.6 > \/2 
[and hence V^tiq) > 17(0), see eq. (|5.1|) ] and g = 10“^. We hnd the typical behaviour of 
unbroken symmetry. Notice again that the ehective or tree level potential is an irrelevant 
quantity for the dynamics, the asymptotic amplitude of oscillation of the zero mode is 
rj ^ 0.5, which is smaller than the minimum of the tree level potential p = 1 but the 
oscillations are symmetric around g = 0. 

Since the dynamical evolution sampled both vacua symmetrically from the beginning, 
there never was a symmetry breaking in the hrst place, and “symmetry restoration” is 
just the statement that the initial state has enough energy density such that the dynamics 
probes both vacua symmetrically despite the fact that the tree level potential allows for 
broken symmetry ground states. 


VII. LINEAR VS. NONLINEAR DISSIPATION (THROUGH PARTICLE 

CREATION) 


As already stressed, the held theory dynamics is unavoidable nonlinear for processes 
like preheating and reheating. It is however interesting to study such processes in the 
amplitude expansion. This is done in detail in refs. |^-[I0| . To dominant order, the amplitude 


expansion means to linearize the zero mode evolution equations. This approach permits an 
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analytic resolution of the evolution in closed form by Laplace transform. Explicit integral 
representations for r]{t) follow as functions of the initial data M- Moreover, the results 
can be clearly described in terms of S-matrix concepts (particle poles, production thresholds, 
resonances, etc.). 

Let us consider the simplest model where the inflaton <h couples to a another scalar a 
and to a fermion held '0, and potential [T^ 


1/ = ^ n7.|<h^ + m^cr^ + g 
In the unbroken symmetry case (m| > 0) the inhaton is always stable and we found for the 


order parameter (expectation value of <h) evolution in the amplitude expansion ||10||, 
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ujT,i{uj) cos cut duj 


lm^+ 2 m„ [a;2 - m| - + T,i{ujy 


(7.1) 


where Sphysicai(*a; ± O’*") = S/j(a;) ± is the inhaton self-energy in the physical sheet, 

Vi = v{^) and v{0) = 0. The hrst term is the contribution of the one-particle pole (at the 
physical inhaton mass). This terms oscillates forever with constant amplitude corresponding 
to the asymptotic single particle state. The second term is the cut contribution v{'t)cut 
corresponding to the process ^ ^ + 2a above the three particle threshold. 

In general, when 

the the cut contribution behaves for late times as 


V{t)cut ^ 
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Here, (Uth = + 2M„ is the threshold energy and a = 2 since to two-loops, 0 




2g^7i^ 


{AtiY (m$ 2M^)^/2 


[uj - (m$ 2MY) ] 


2i2 


In the broken symmetry case (m| < 0) we may have either M < 2mo- or M > 
where M is the physical inhaton mass. ( M = |m$|A/2 at the tree level). In the hrst case 


the inhaton is stable and eq. (|7.1|) holds. However, the self-energy starts now at one-loop 
and vanishes at threshold with a power a = 1/2. For M > the “inhaton” becomes a 
resonance (an unstable particle) with width (inverse lifetime) 


T = 
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This pole dominates r]{t) for non asymptotic times 


5{t) ^ 6iA e cos(Mt + 7) , 


(7.2) 


where 


71 = 1 + 




7 = 


dEi{M) 


dM^ ’ ' dM^ ■ 

In summary, eq. (|7.2|) holds provided: a) the inflaton is a resonance and b) t < T”^ ln(r/Mo-). 
For later times the fall off is with a power law determined by the spectral density at 
threshold as before |T^. The full study of “inflaton” decaying into lighter scalars in DeSiter 
space-time has been recently presented in reference 


In ref. the selfconsistent nonlinear evolution is computed to one-loop level for the 
model (|7.1| ). In £g. 8 ri{t) is plotted as a function of time ior X = g = l.bvr^, y = 0, m„ = 
0 . 2 m$, 7 ( 0 ) = 1 and 7 ( 0 ) = 0 . 

Figure 8 shows a very rapid, non-exponential damping within few oscillations of the 
expectation value and a saturation effect when the amplitude of the oscillation is rather 
small (about 0.1 in this case), the amplitude remains almost constant at the latest times 
tested. Figures 8 and 9 clearly show that the time scale for dissipation (from £g. 8 ) is 
that for which the particle production mechanism is more efficient (fig. 9). Notice that the 
total number of particles produced rises on the same time scale as that of damping in fig. 8 
and eventually when the expectation value oscillates with (almost) constant amplitude the 
average number of particles produced remains constant. This behaviour is a close analog to 
the selfcoupled inflaton for unbroken symmetry (fig.l). The amplitude expansion predictions 
are in qualitative agreement with both results. 

These figures clearly show that damping is a consequence of particle production. At times 
larger than about 40 (for the initial values and couplings chosen) there is no appreciable 
damping. The amplitude is rather small and particle production has practically shut off. If 
we had used the classical evolution of the expectation value in the mode equations, particle 
production would not shut off (parametric resonant amplification), and thus we clearly see 
the dramatic effects of the inclusion of the back reaction. 

In ref. the broken symmetry case m|, < 0 is then studied. Figures 11-13 show t]{t) 
vs r, Mair) vs r and = 200) vs q respectively, for X/Sn"^ = 0 . 2 ; g/X = 0.05; = 

0.2 |m$|; 7 ( 0 ) = 0.6; r]{0) = 0. Notice that the mass for the linearized perturbations of the 
$ field at the broken symmetry ground state is ^/2\mq>\ > 2m„. Therefore, for the values 
used in the numerical analysis, the two-particle decay channel is open. For these values of the 
parameters, linear relaxation predicts exponential decay with a time scale Trei ~ 300 (in the 
units used). Figure 11 shows very rapid non-exponential damping on time scales about six 
times shorter than that predicted by linear relaxation. The expectation value reaches very 
rapidly a small amplitude regime, once this happens its amplitude relaxes very slowly. In the 
non-linear regime relaxation is clearly not exponential but extremely fast. The amplitude 
at long times seems to relax to the expected value, shifted slightly from the minimum of 
the tree level potential at 7 = 1. This is as expected from the fact that there are quantum 
corrections. Figure 12 shows that particle production occurs during the time scale for which 
dissipation is most effective, giving direct proof that dissipation is a consequence of particle 
production. Asymptotically, when the amplitude of the expectation value is small, particle 
production shuts off. We point out again that this is a consequence of the back-reaction in 
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the evolution equations. Without this back-reaction, as argued above, particle production 
would continue without indehnitely. Figure 13 shows that the distribution of produced 
particles is very far from thermal and concentrated at low momentum modes k < |m$|. 
This distribution is qualitatively similar to that in the unbroken symmetry case, and points 
out that the excited state obtained asymptotically is far from thermal. 

In ref. 13 the case where the inflaton is only coupled to fermions is studied {g = 0, y ^ 
0). The damping of the zero mode is very inefficient in such case due to Pauli blocking. 
Namely, the Pauli exclusion principle forbids the creation of more than 2 fermions per 
momentum state. Pauli blocking shuts off particle production and dissipation very early on. 


VIII. THE REHEATING TEMPERATURE 


The arena in which these results become important is that of inflationary cosmology. 
In particular, the process of preheating is of vital importance in understanding how the 
big-bang cosmology is regained at the end of inflation, i.e. the reheating mechanism. 

While our analysis has been entirely a Minkowski space one, we can make some comments 
concerning the reheating temperature. However, a more detailed analysis incorporating the 
expanding universe must eventually be done along the lines suggested in this work, to get 
more accurate results. 

Since the particles created during the preheating stage are far from equilibrium, thermal- 
ization and equilibration will be achieved via collisional relaxation. In the approximation 
that we are studying, however, collisions are absent and the corresponding contributions are 
of 0(1/N). The difficulty with the next order calculation and incorporation of scattering 
terms is that these are non-local in time and very difficult to implement numerically. 

However we can obtain an estimate for the reheating temperature under some reasonable 
assumptions: in the cosmological scenario, if the equilibration time is shorter than the inverse 
of the expansion rate H, then there will not be appreciable redshifting of the temperature 
because of the expansion and we can use our Minkowski space results. 

The second assumption is that the time scales between particle production and ther- 
malization and equilibration are well separated. Within the large N approximation this is 
clearly correct because at large N, scattering processes are suppressed by 1/N. If these 
two time scales are widely separated then we can provide a fairly reliable estimate of the 
reheating temperature as follows. 

Equilibration occurs via the redistribution of energy and momentum via elastic collisional 
processes. Assuming that thermalization occurs on time scales larger than that of particle 
production and parametric amplihcation, then we can assume energy conservation in the 
scattering processes. Although a reliable and quantitative estimate of the reheating temper¬ 
ature can only be obtained after a detailed study of the collisional processes which depend 
on the interactions, we can provide estimates in two important cases. If the scattering pro¬ 
cesses do not change chemical equilibrium, that is conserve particle number, the energy per 
particle is conserved. Since the energy (density) stored in the non-equilibrium bath and the 
total number of particles per unit volume are, as shown in the previous sections: 


\Mr 

A 
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N 







with proportionality constants of order one, we can estimate the reheating temperatnre to 
be: 

±^T^\Mn\ ( 8 . 1 ) 

Here \Mji\ is the inflaton mass. This is consistent with previous results Q. 

This result seems puzzling, since naively one would expect £ ~ T^; N ^ but the 
powers of A do not match. This puzzle arises from intuition based on an ultrarelativistic 
free particle gas. However the “medium” is highly excited with a large density of particles 
and the “in medium” properties of the equilibrated particles may drastically modify this nave 
result as is known to happen in most theories at high temperature, where the medium effects 
are strong and perturbation theory breaks down requiring hard thermal loop resummation. 

In the case in which the collisional processes do not conserve particle number and there¬ 
fore change chemical equilibrium, the only conserved quantity is the energy. Such is the 
case for massless particles interacting with a quartic couplings for example or higher order 
processes in a quartic theory with massive particles. Processes in which 3 —> 1 conserving 
energy and momentum can occur. The inverse process 1 —> 3 occurs with far less probability 
since the high momentum modes are much less populated than the low momentum modes 
in the unstable bands. In this case only energy is conserved whereas the total number of 
particles is not conserved and in this case an estimate of the reheating temperature compares 
the energy density in the bath of produced particles to that of an ultrarelativistic gas in 
equilibrium at temperature Tr, 


leading to the estimate 


\Mr[- 


A 


R 


-R 


Tr 


\Mr\ 

Ar-* 


Thus we can at least provide a bound for the reheating temperature 


|aPr| <Tr< 


\Mr\ 

A 


and a more quantitative estimate requires a deeper understanding of the collisional processes 
involved. 

Within the large N approximation, scattering terms will appear at order 1/N and beyond. 
The leading contribution 0{1/N) to collisional relaxation conserves particle number in the 
unbroken symmetry state because the product particles are massive. This can be seen from 
the fact that the self-energy to this order is given by the same chain of bubbles that gives 
the scattering amplitude but with two external legs contracted. All cut diagrams (that 
give the imaginary part) correspond to 2 —>■ 2 processes that conserve particle number 
because of kinematic reasons. Certainly at higher order in 1/A^ there will be processes that 
change chemical equilibrium, but for large N these are suppressed formally. This argument 
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based on the leading collisional contribution in the 1/N expansion allows us to provide a 
further consistent estimate in the unbroken symmetry case (when the produced particles 
are massive). In this approximation and consistently with energy and particle number 
conservation we can assume that the hnal equilibrium temperature is of the order of the 
typical particle energy before thermalization. Recalling that the unstable band remains 
stable during the evolution with the peak shifting slightly in position we can estimate the 
typical energy per particle by the position of the peak in the distribution at gi, and use the 
analytical estimate for the peak given in section 2. Restoring the units we then obtain the 
estimate 


Tr k. \Mr\ qi 





( 8 . 2 ) 


which displays the dependence on rjo explicitly. Eq. ( p.2|) is an improvement of the simple 
estimate S). 

It must be noticed that the peak in the momentum distribution decreases with time for 
t > treh (see £g.2f-2h). This drift follows from the non-linear interaction between the modes. 
For the case of £g. 2, one sees that Tr reduces by approximately a factor 3 with respect to 
the value (|8.2|) . 

The large N model studied in this article is not a typical model used in inflationary 
cosmology, and since we want to make a quantitative statement for inflationary scenarios 
(within the approximation of only considering Minkowski space) we now study a model that 
incorporates other scalar helds coupled to the inflaton. 

The simplest model [@] contains in addition to the inflaton a lighter scalar held a with a 
coupling. That is, we consider the Lagrangian 0], 


C = “ 2 *^ + rri^ + ga) ^ 


4! 


$4 


^ /o2 I 2 \ 4 

-a(d +m^)a-—a 


We will consider again the preheating regime of weak couplings and early times such 
that we can neglect the back-reaction of the quantum huctuations of the a held as well 
as the back-reaction of the quantum huctuations of the inhation itself, focussing only on 
the parametric growth of the a huctuations in the unbroken symmetry case. The mode 
equations for the a held take the form 




In dimensionless variables this equation becomes 


fk(t) = 0 . 




/m„ 


^ + q +{— +^v{r 




\ m 


6g 


X 


fkir) = 0 . 


In the short time approximation we can replace //(r) by the classical form (|4.1| ). We then 
hnd a Lame equation which admits closed form solutions for 


A 


= n{n + 1) 


n = 1, 2, 3,... 
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Although these are not generic values of the couplings, the solubility of the model and 
the possibility of analytic solution for these cases makes this study worthwhile. In the 
simplest case, (n = 1, 6g = X), there is only one forbidden band for > 0. It goes from 

q^ = l-(^^y to q^ = l-{^y + iql/2. That is. 


forbidden band : m? — mt < 


4 




The Floquet index in this forbidden band takes the form 

F{q) = 2i K{k) Z{2K{k) n) ± tt , 


where now v and q are related by the following equation 

\ m J 2 

and 0 < n < 1/2. 

The imaginary part of the Floquet index is now maximal at = 1 — ^ and we 

can use this value to provide an estimate for the reheating temperature in this model in the 
same manner as for eq.(^.2|), yielding the following estimate for the reheating temperature 


Treh ^ 




rrii + 




(8.3) 


In the special case = m = we recover eq.(^.2D, as expected. 

Again, the non-lineal field evolution for t > treh decreases Treh- In the third reference 
under [^, we found for late times a Treh ten times smaller than the value ( |8.3|) for g = l.bvr^. 

The estimates on the reheating temperature provided above should not be taken rigor¬ 
ously, but as an approximate guide. A consistent estimate of the reheating temperature and 
the thermalization time scales would in principle involve setting up a Boltzmann equation 
0]. Under the assumption of a separation between the preheating and thermalization time 
scales one could try to use the distribution functions Ng at the end of the preheating stage 
as input in the kinetic Boltzmann equation. However, we now argue that such a kinetic de¬ 
scription is not valid to study thermalization. A kinetic approach based on the Boltzmann 
equation, with binary collisions, for example, would begin by writing the rate equation for 
the distribution of particles 


iVfc oc A^y d^kid?k 2 d^k-i 5^{ki + k 2 + k-i + k) [(1 nk){l + nfcJnfcjUfcg - 

nknk^{l + nk^){l + Uk^)] . 


However this equation is only valid in the low density regime. In particular for the case under 
study, the occupation numbers for wavevectors in the unstable bands are non-perturbatively 
large oc 1/A and one would be erroneously led to conclude that thermalization occurs on the 
same time scale or faster than preheating. 

Clearly such a statement would be too premature. Without a separation of time scales 
the kinetic approach is unwarranted. The solution of the Boltzmann equation provides a 
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partial resummation of the perturbative series which is valid whenever the time scales for 
relaxation is much longer than the microscopic time scales [^|, in this case that of particle 
production. 

In the case under study there is an expansion parameter, 1/N and clearly these scat¬ 
tering terms are sub leading in this formal limit, so that the separation of time scales is 
controlled. In the absence of such an expansion parameter, some resummation scheme must 
be invoked to correctly incorporate scattering. In particular when the symmetry is broken, 
the asymptotic excitations are Goldstone bosons, the medium is highly excited but with 
very long-wavelength Goldstones and these have very small scattering cross sections. Such a 
resummation is also necessary in the large temperature limit of held theories in equilibrium. 
In this case the perturbative expansion of the scattering cross section involves powers of 
XT/m with m being the mass. A correct resummation of the (infrared) divergent terms 
leads to A(T) —>■ m/T in the large T/m limit 


37|| . In particular the 1/N corrections in the 


formal large N limit involve such a resummation, but in the non-equilibrium situation, the 
numerical implementation of this resummation remains a formidable problem. 


IX. CONCLUSIONS 

It is clear that preheating is both an extremely important process in a variety of settings, 
as well as one involving very delicate analysis. In particular, its non-perturbative nature 
renders any treatment that does not take into account effects such as the quantum back- 
reaction due to the produced particles, consistent conservation (or covariant conservation) of 
the relevant quantities and Ward identities, incapable of correctly describing the important 
physical phenomena during the preheating stage. 

In this work, we dealt with these issues by using the 0{N) vector model in the large N 
limit. This allows for a controlled non-perturbative approximation scheme that conserves 
energy and the proper Ward identities, to study the non-equilibrium dynamics of scalar 
helds. Using this model we were able to perform a full analysis of the evolution of the zero 
mode as well as of the particle production during this evolution. 

Our results are rather striking. We were able to provide analytic results for the held 
evolution as well as the particle production and the equation of state for all these components 
in the weak coupling regime and for times for which the quantum huctuations, which account 
for back-reaction effects, are small. What we found is that, in the unbroken symmetry 
situation, the held modes satisfy a Lame equation that corresponds to a Schrodinger equation 
with a two-zone potential. There are two allowed and two forbidden bands, which is decidedly 
unlike the Mathieu equation used in previous analysis The diherence between an 

equation with two forbidden bands and one with an inhnite number is profound. We were 
also able to estimate analytically the time scale at which preheating would occur by asking 
when the quantum huctuations as calculated in the absence of back-reaction would become 
comparable to the tree level terms in the equations of motion. The equations of state of both 
the zero mode and “pions” were calculated and were found to be describable as polytropes. 
These results were then conhrmed by numerical integration of the equations, and we found 
that the analytic results were in great agreement with the numerical ones in their common 
domain of validity. 
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When the 0{N) symmetry is spontaneously broken, more subtle effects can arise, again 
in the weak coupling regime. If the zero mode starts off very near the origin, then the 
quantum back-reaction grows to be comparable to the tree level terms within one or at most 
a few oscillations even for very weak coupling. In this case the periodic approximation for 
the dynamics of the zero mode breaks down very early on and the full dynamics must be 
studied numerically. 

When numerical tools are brought to bear on this case we hnd some extremely interesting 
behavior. In particular, there are situations in which the zero mode starts near the origin 
(the initial value depends on the coupling) and then in one oscillation, comes back to almost 
the same location. However, during the evolution it has produced 1/g particles. Given that 
the total energy is conserved, the puzzle is to hnd where the energy came from to produce the 
particles. We found the answer in a term in the energy density that has the interpretation 
of a “vacuuum energy” that becomes negative during the evolution of the zero mode and 
whose contribution to the equation of state is that of “vacuum”. The energy given up by 
this term is the energy used to produce the particles. 

This example also allows us to study the possibility of symmetry restoration during 
preheating While there have been arguments to the effect that the produced 

particles will contribute to the quantum fluctuations in such a way as to make the effective 
mass squared of the modes positive and thus restore the symmetry, we argued that they 
were unfounded. Whereas the effective squared mass oscillates, taking positive values during 
the early stages of the evolution, its asymptotic value is zero, compatible with Goldstone 
bosons as the asymptotic states. 

Furthermore, this says nothing about whether the symmetry is restored or not. This is 
signaled by the hnal value of the zero mode. In all the situations examined here, the zero 
mode is driven to a non-zero hnal value. At this late time, the “pions” become massless, i.e. 
they truly are the Goldstone modes required by Goldstone’s theorem. 

The arguments presented in favor of symmetry restoration rely heavily on the effective 
potential. We have made the point of showing explicitly why such a concept is completely 
irrelevant for the non-equilibrium dynamics when profuse particle production occurs and 
the evolution occurs in a highly excited, out of equilibrium state. 

Finally, we dealt with the issue of how to use our results to calculate the reheating 
temperature due to preheating in an inhationary universe scenario. Since our results are 
particular to Minkowski space, we need to assume that preheating and thermalization occur 
on time scales shorter than the expansion time, i.e. H~^. We also need to assume that 
there is a separation of time scale between preheating and thermalization. Under these 
assumptions we can estimate the reheating temperature as T^eh oc \Mp\ in the case where 
the produced particles are massive and T^eh c< |Mj^|/A4 in the massless case. 

We have made the important observation that due to the large number of long-wavelength 
particles in the forbidden bands, a kinetic or Boltzmann equation approach to thermalization 
is inconsistent here. A resummation akin to that of hard thermal loops, that consistently 
arises in the next order in 1/A^ must be employed. In equilibrium such a resummation shows 
that the scattering cross section for soft modes is perturbatively small despite their large 
occupation numbers. 

There is a great deal left to explore. Recently we reported on our study of the non¬ 
linear quantum held evolution in de Sitter and FRW backgrounds in refs. and ||17||, 
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respectively. The next important step is to consider the background dynamics as a full 
backreaction problem, including the inflaton dynamics and the dynamics of the scale factor 
self-consistently. 

Such a detailed study will lead to a consistent and thorough understanding of the infla¬ 
tionary period, the post-inflationary and reheating periods. Further steps should certainly 
include trying to incorporate thermalization effects systematically within the 1/N expansion. 


The preheating and reheating theory in inflationary cosmology is currently a very active 
area of research in fast development, with the potential for dramatically modifying the 
picture of the late stages of inflationary phase transitions. 

As remarked before, reliable and consistent estimates and held theory calculations have 
been done mostly assuming Minkowski spacetime. The matter state equations obtained 
in Minkowski and recently in de Sitter backgrounds |]^ give an indication, through 
the Einstein-Friedmann equation, of the dynamics of scale factor and give a glimpse of the 
important physics to be unraveled by a deeper study. 

The formulation described in detail in these lectures are uniquely suited to provide 
complete description of the full dynamics of inhationary cosmology, from times prior to 
the phase transitions or the beginning of the chaotic era, through the inhationary regime, 
to the post-inhat ionary and reheating stage. 

Such a program provides the ultimate tool to test physical predictions of particle physics 
models. Thus this new consistent formulation provides the practical means to input a 
particle physics model and extract from it reliable dynamical predictions which will have to 
ultimately be tested against the next generation of cosmological experiments. 
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TABLES 


Vo 

q 

B 

N 

1 

0.017972387... 

0.1887167... 

3.778... 

3 

0.037295557... 

0.8027561... 

0.623... 

4 

0.03966577... 

1.1007794... 

0.4... 

r]o^ oo 

0.043213918... 

0.28595318 r/o + 0(%^) 

3.147r/o'/'[l + 0(r,o-2)] 


TABLE 1. Quantum Fluctuations S(r) ]v^ during the preheating period. 
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Figure 1: The ratio < po > /^o for zero mode vs. XReo/2\Mfi\‘^ for the unbroken 
symmetry case. 
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Figure 2 b 



Figure 2b: r]{T) for the same values of the parameters as in Fig. 2(a). The agreement 
with the analytic prediction is to within 5% for 0 < r < 30. 
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Figure 2c: ( 7 S(r) for the same parameters as in fig. 2a. 














































































































































































































































Figure 2d: g'E{T), analytic approximation and nnmerical resnlts for the same valnes of 
parameters as in fig. 2a. 


Figure 2 d 





Figure 2e: gNij) vs. r for the same values of parameters as fig. 2(a). 






























Figure 2f: gNq{T) vs. q for r = 200 for the same values of parameters as fig. 


Figure 2 f 
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totically the average over a period gives p 


















































































































Figure 4b: 5 'S(r) for the same values of the parameters as in hg. 4a. 










Figure 4c: gM{T) for the same parameters as in fig. 4a. 
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Figure 4d: gNqij) vs. q for r = 30 for the same values of parameters as Fig. 4a. 


Figure 4 d 








Figure 4e: gNq{T) vs. g for r = 90 for the same values of parameters as Fig. 4a. 


Figure 4 e 








Figure 4f: gNq{T) vs. q for r = 150 for the same values of parameters as Fig. 4a. 
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Figure 4g: A^^(r) vs. r for the same parameters as Fig. 4a. 
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Figure 4i: eiv(r) vs. r for the same parameters as Fig. 4a. 











Figure 4j: £c(t) vs. r for the same parameters as Fig. 4a. 
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Figure 5b: A^^(r) vs. r for the same parameters as fig. 5(a). 
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Figure 7: The effective mass squared vs. r for the same values of the parameters as in 






















87 









































90 








